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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
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Tennessee.  The  work  was  done  under  ARO  Project  Nos.  P33A-36A  and  P32A-31D.  The 
author  of  this  report  was  David  L.  Whitfield,  ARO,  Inc.  The  manuscript  (ARO  Control 
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1.0  INTRODUCTION 

The  material  presented  herein  is  a  result  of  work  done  in  support  of  various  research, 
test,  and  test  facility  development  projects  conducted  in  the  Propulsion  Wind  Tunnel 
Facility  (PWT).  The  topics  considered  are  somewhat  independent,  and  an  effort  is  made 
to  develop  each  in  an  essentially  self-contained  section.  This  permits  the  results  to  be 
reviewed  and  extracted  for  particular  applications  as  desired  without  sorting  through  the 
entire  report.  However,  an  important  use  of  these  results  is  their  application  as  a  unit 
to  the  analysis  of  experimental  data.  Such  an  application  is  made,  not  only  as  an  example 
to  illustrate  this  point,  but  to  analyze  experimental  turbulent  boundary-layer  measurements 
made  in  PWT. 

The  first  topic  considered  (Section  2.0)  is  a  laminar  and/or  turbulent  numerical 
boundary-layer  calculation  technique  in  which  the  Reynolds  stress  is  related  to  the  local 
turbulent  kinetic  energy  of  the  flow,  e,  and  a  modified  form  of  the  turbulent  kinetic 
energy  equation  is  used  to  solve  for  e  throughout  the  boundary  layer.  This  approach  has 
certain  advantages  over  the  classical  mixing-length  approach.  For  example,  computations 
can  be  performed  of  transition  flow,  initial  and  free-stream  turbulence  levels  can  be  taken 
into  account,  and  the  integration  is  carried  completely  to  the  wall  where  natural  boundary 
conditions  can  be  imposed.  This  numerical  calculation  technique  is  developed  for  the 
purpose  of  performing  reasonable  engineering  predictions  of  incompressible  and 
compressible  turbulent  boundary  layers.  The  major  effort  of  this  portion  of  the  report 
is  in  the  modeling  of  the  various  terms  of  the  turbulent  kinetic  energy  equation. 

Section  3.0  contains  the  development  of  an  analytical  investigation  of  turbulence 
near  a  wall.  The  approach  is  unique  in  that  it  is  based  on  the  turbulent  kinetic  energy 
equation  developed  in  Section  2.0,  and  it  does  not  depend  on  mixing-length  or 
damping-factor  concepts.  The  results  include  closed-form  expressions  for  the  velocity, 
Reynolds  stress,  production  of  turbulence,  and  dissipation  of  turbulence.  Moreover,  the 
results  show  good  agreement  with  experimental  data  and  are  valid  from  the  wall  through 
the  so-called  sublayer  and  buffer  layer  and  into  the  fully  turbulent  portion  of  the  boundary 
layer.  An  important  application  of  these  results  is  that  of  providing  information  of  the 
region  near  the  wall  where  experimental  measurements  are  not  routinely  made  (usually 
never  made  during  wind-tunnel  tests  because  of  the  experimental  difficulties  and/or  data 
acquisition  time).  These  analytical  results  are  used  to  more  accurately  determine 
boundary-layer  parameters  from  experimental  data  such  as  displacement  thickness  and 
momentum  thickness  as  discussed  in  Section  5.0. 
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Section  4.0  is  concerned  with  the  development  of  analytical  relationships  between 
velocity  and  temperature  throughout  turbulent  boundary  layers  on  constant  temperature 
or  adiabatic  walls  for  nonunity  Prandtl  numbers.  Analytical  solutions  are  of  interest  as 
opposed  to  numerical  ones  for  purposes  of  data  reduction  of,  e.g.  pitot  pressure 
measurements  and  holographic  interferometric  measurements.  In  addition,  Reynolds 
analogy  is  calculated  from  these  solutions  which  is  shown  to  be  in  good  agreement  with 
previous  results  for  near  adiabatic  wall  conditions,  and  predicts  the  quantity  2  Ch/cf  to 
decrease  with  decreasing  wall  temperature.  This  prediction  is  also  shown  to  be  supported 
by  comparisons  with  experimental  data,  and  apparently  such  a  variation  has  not  previously 
been  predicted. 

In  Section  5.0  the  result  for  the  velocity  distribution  and  the  velocity-temperature 
relations  (Sections  3.0  and  4.0)  are  included  in  the  development  of  a  data  reduction 
computer  program  for  the  purpose  of  reducing  pitot  pressure  measurements  made  in 
turbulent  boundary  layers.  The  code  is  written  to  obtain  spatial  distributions  of  velocity 
and  temperature,  local  skin  friction,  and  various  boundary-layer  parameters. 

The  results  of  Sections  2.0  through  5.0  are  applied  in  Section  6.0  to  the  analysis 
of  turbulent  boundary-layer  measurements  made  in  Propulsion  Wind  Tunnels  16S,  16T, 
and  the  Acoustic  Research  Tunnel  (ART)  in  PWT.  The  experimental  data  considered  include 
previously  reported  measurements  made  in  Tunnel  1 6S  and  recently  acquired  unpublished 
data  in  Tunnels  16S,  16T,  and  ART.  Results  of  this  analysis  which  are  of  particular  interest 
include  (1)  new  skin-friction  coefficients  determined  from  measurements  made  in  Tunnel 
16S  and  their  impact  on  transition  correlations  of  data  taken  in  16S  and  (2)  the  nature 
of  the  flow  in  the  nozzle  boundary  layer  of  Tunnel  16T. 

This  report  contains  a  fair  amount  of  new  results.  For  example,  (1)  the  particular 
modeling  of  the  turbulent  kinetic  energy  equation  developed  in  Section  2.0  has  not  been 
used  previously,  (2)  the  analytical  results  obtained  in  Section  3.0  concerning  turbulence 
near  a  wall  are  new,  but  more  importantly  this  particular  development  (i.e.  without 
mixing-length  theory)  has  not  been  followed  previously,  and  (3)  the  analytical 
velocity-temperature  relations  and  Reynolds  analogy  for  nonunity  Prandtl  numbers 
obtained  in  Section  4.0  are  also  new.  Therefore,  because  these  results  are  not  common 
nor  time  tested,  effort  was  expended  in  assessing  the  accuracy  of  each  by  numerous 
comparisons  with  experimental  data.  In  Sections  2.0  through  4.0,  these  experimental 
comparisons  follow  the  theoretical  development  of  the  results.  The  computer  programs 
for  the  boundary-layer  calculations  and  the  data  reduction  procedure  (Sections  2.0  and 
5.0)  are  included  in  Appendixes  A  and  B,  respectively. 
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2.0  NUMERICAL  BOUNDARY-LAYER  CALCULATIONS 

This  numerical  technique  draws  heavily  from  that  of  previous  investigators.  Portions 
of  several  sources  have  been  extracted  and  modified  where  deemed  necessary  for  the 
purpose  of  improving  the  results.  The  objective  was  to  obtain  as  accurate  results  as  possible 
(the  standard  for  which  was  comparisons  with  experimental  data)  for  various  flow 
conditions,  while  attaining  a  physically  reasonable  model  with  minimum  empiricism  (the 
constants  of  which  were  not  allowed  to  change  for  the  various  flow  conditions).  The 
result  is  a  fairly  easy  to  use  numerical  method  which  yields  results  in  relatively  good 
agreement  with  experimental  data  for  the  flows  considered.  However,  there  is  much  room 
for  improvement.  For  example,  the  present  code  will  run  for  laminar  or  turbulent  flow 
and  calculate  through  transition,  producing  quite  reasonable  looking  velocity  profiles. 
Unfortunately,  the  point  of  transition  occurs  upstream  (lower  Reynolds  number)  of  that 
observed  experimentally.  Little  has  been  attempted  with  the  present  code  concerning 
possible  improvements  in  calculating  transition. 

It  should  be  noted  that  this  investigation  was  made  much  easier  by  the  powerful 
numerical  method  for  diffusion-type  equations  developed  by  Patankar  and  Spalding  (Ref. 
1).  Effort  here  was  expended  in  modeling  the  equations  rather  than  in  their  solution. 
A  listing  of  the  resulting  computer  code  is  given  in  Appendix  A. 

2.1  DIMENSIONLESS  BOUNDARY-LAYER  EQUATIONS 

The  choice  of  the  dimensionless  variables  was  such  as  to  recover  the  general  form 
of  the  physical-plane  equations  considered  by  Patankar  and  Spalding  (Ref.  1)  who  apphed 
the  von  Mises  transformation  to  the  physical-plane  equations,  whereas  here  the  von  Mises 
transformation  was  applied  to  dimensionless  equations.  Dimensionless  variables  simplify 
and  reduce  appreciably  the  input  required  to  run  the  program.  Also  some  physical  insight 
into  the  flow  problem  and  its  relation  to  other  flows  is  gained  by  using  dimensionless 
parameters  and  eliminating  the  confusion  attributable  to  units. 

In  the  dimensionless  variables  defined  in  the  nomenclature  and  used  in  Ref.  2,  the 
boundary-layer  equations  considered  were 

dip  r  U)  ^  d(p  r  V)  ^  ^ 

pu 


d\]  -tt 

—  +  p  V  —  = 

dx  dy 


o,  s 


_  ^  \ 

dU/dy  dy  | 


p  <uv> 


(2) 
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r  =  ±  y  cos  0 


(5) 


where  the  plus  sign  in  Eq.  (5)  is  for  external  flow  and  the  minus  sign  for  internal  flow. 
The  quantity  (2)1  /Kqq  as  a  coefficient  of /i  in  Eqs.  (2)  and  (3)  is  due  to  using  (Hq 
to  nondimensionalize  velocities  rather  than  (2Ho,^)i/^,  the  maximum  velocity.  Using 
(2Ho,^)i/^  would  introduce  new  coefficients  for  other  terms,  which  were  thought  to  be 
less  convenient  to  handle.  The  (2)1 was  used  in  Reo,s  simply  for  convenience  in  order 
that  Reo,s  based  on  the  maximum  velocity. 


In  order  to  close  this  system  of  equations  something  must  be  done  with  the  Reynolds 
stress  term.  Rather  than  introduce  the  usual  eddy-viscosity  or  mixing-length  concepts,  the 
Reynolds  stress  was  modeled  here  by  relating  it  to  the  turbulent  kinetic  energy. 


2.2  MODEL  FOR  REYNOLDS  STRESS 


Two  forms  are  commonly  used  for  modeling  the  Reynolds  stress  by  relating  it  to 
the  turbulent  kinetic  energy  of  the  flow.  The  first,  introduced  by  Townsend  (Ref.  3) 
and  refined  by  Lighthill  (Ref.  4),  written  in  the  boundary-layer  approximation,  is 

-  <\iv>  =  ae  (6) 

and  the  second,  introduced  by  Prandtl  (Ref.  5),  is 

-  <uv>  =  (7) 

dy 

where  a  and  b  are  constants  and  £  is  the  mixing  length.  Equation  (6),  with  a  =  0.3, 
was  used  in  the  present  work.  Both  equations  have  been  used  by  previous  investigators. 
For  example,  Bradshaw,  Ferriss  and  Atwell  (Ref.  6)  used  Eq.  (6)  with  a  =  0.3;  Ng  and 
Spalding  (Ref.  7)  used  Eq.  (6)  in  equilibrium  flow  for  evaluating  a  constant  in  the 
dissipation  term  (a  effectively  was  0.32);  and  Glushko  (Ref.  8)  (recently  extended  by 
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Beckwith  and  Bushnell  (Ref.  9))  used  a  model  which  can  be  reduced  to  Eq.  (7)  over 
a  portion  of  the  boundary  layer.  These  approaches,  as  well  as  the  present  one,  differ 
primarily  in  the  modehng  of  the  terms  in  the  turbulent  kinetic  energy  equation  which 
is  required  for  a  solution  for  e. 

2,3  TURBULENT  KINETIC  ENERGY  EQUATION 


The  derivation  of  the  particular  turbulent  kinetic  energy  equation  used  here  begins 
with  the  equation  proposed  by  Kolmogoroff  (Ref.  10),  and  used  by  Ng  and  Spalding 
(Ref.  7),  and  is  written  for  axisymmetric  flow  as 


'\T  ^  ^ 


J_  I  ± 

(7g  r 


(8) 


The  terms  in  Eq.  (8)  on  the  left-hand  side  represent  convection,  and  the  first,  second, 
and  third  terms  on  the  right-hand  side  represent  the  diffusion,  production,  and  dissipation, 
respectively.  Glushko  (Ref.  8)  points  out  that  at  low  Reynolds  numbers  the  turbulent 
energy  dissipation  deviates  from  that  specified  in  Eq.  (8).  The  reason  for  this  has  to  do 
with  assumptions  concerning  the  equilibrium  of  the  turbulence  at  high  wave  numbers  for 
sufficiently  high  Reynolds  numbers  as  explained  by  Batchelor  (Ref.  11).  Ng  and  Spalding 
(Ref.  7)  used  Eq.  (8)  in  conjunction  with  an  additional  equation  in  order  to  calculate 
However,  besides  the  problem  with  the  dissipation  term,  Ng  and  Spalding  (Ref.  7) 
did  not  use  the  molecular  viscosity  and  did  not  integrate  completely  to  the  wall.  In  an 
attempt  to  remedy  these  deficiencies,  the  following  modifications  were  made  to  Eq.  (8). 
The  numerator  and  denominator  of  the  dissipation  term  in  Eq.  (8)  were  multiplied  by 
the  integral  scale  length  of  turbulence  £',  in  order  to  obtain  a  term  e^/^  £'  common 
to  each  term  on  the  right-hand  side  of  Eq.  (8).  Ng  and  Spalding  (Ref.  7)  used  a  proposal 
of  Kolmogoroff  (Ref.  10)  and  Prandtl  (Ref.  5)  to  relate  e^/^  £'  to  the  Reynolds  stress 
by  el/2  Q'  =  _  <uv>/3U/9y  which  is  the  form  of  Eq.  (7).  However,  in  the  present  work 
the  expression 


e 


^2 


<UV>  fl 

idU/dy)  ^  p 


(9) 


is  used  in  the  diffusion  and  dissipation  terms  while  retaining  ei/2  £'  =  -  <uv>/dU/3y 
in  the  production  term.  Equation  (8)  can  then  be  written 


de  de  1 


d  /  p<nv>\  de 
dy  ^  \  (9U /dy  j  dy 


.  ^  (9U 
p<iiv>  -Tj— 
dy 


Ce  /  p<uv>\ 
^2  d\]/dyj 


(10) 


where  L  is  the  dissipation  length  and  C  is  a  constant  or  a  function  of  Reynolds  number. 
Equation  (9)  was  not  used  in  the  production  term  in  Eq.  (8)  because  this  term  represents 
the  production  of  turbulent  energy  and  would  not  be  zero  for  e  =  0  (laminar  flow). 
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Equation  (10)  is  the  turbulent  kinetic  energy  equation  used  here.  It  is  similar  in 
form  to  the  equation  developed  by  Glushko  (Ref.  8).  However,  Glushko  modeled  the 
Reynolds  stress  more  nearly  like  that  given  by  Eq.  (7),  whereas  in  the  present  work  the 
Reynolds  stress- is  modeled  by  Eq.  (6).  Also,  Glushko’s  modeling  takes  on  three  different 
exphcit  forms  across  the  boundary  layer,  whereas  here  Eq.  (7)  was  used  throughout  the 
boundary  layer.  This  modeling  was  used  for  simplicity  and  can  be  justified  only  by  the 
results  obtained.  In  dimensionless  variables  the  additional  equations  used  to  close  the  system 
composed  of  Eqs.  (1)  through  (5)  are 


de  _Yr  de  11^ 

pU—  +pV  —  =  —  z  — 
dx  dy  r  dy 

-  dV 

-  p  <uv>  -3 — 
dy 


(2^  p<5v> 

» 

^®o,s  <9U/^y 

dyj 

r 

“1 

Ce 

l2 


O,  S 


p<uv> 


dM/dy 


(11) 


and 


-  <uv>  =  ae  =  0.3  e  (12) 

The  expression  used  for  C  in  Eq.  (11)  is 

C  =  3.1  -  1.82  X  10-9  (13) 

with  a  lower  bound  of  2.7.  The  reason  for  this  particular  expression  will  be  made  clear 
subsequently  when  comparisons  with  experimental  data  are  discussed. 


Application  of  the  von  Mises  transformation  to  the  boundary-layer  equations  is  given 
in  general  form  by  Patankar  and  Spalding  (Ref.  1)  and  in  the  present  dimensionless 
variables,  except  for  the  additional  terms  due  to  the  Reynolds  stress,  in  Ref.  2.  Applying 
this  transformation  to  the  turbulent  kinetic  energy  equation,  Eq.  (11),  gives  (the  necessary 
partial  derivative  operators  are  given  in  Ref.  2) 


*E  -  n  +  'I  *I  (9e  1  d  \ 

t  ---2 

1  pUr 

(2)"^//  p<uv> 

dt\ 

^  +  4>^  -  da  ~  dcoj 

[(0E  “ 

“'o.s  dv/dy 

dcoi 

(14) 

\?  L^®o,s  dU/dyJ 

The  boundary  conditions  for  the  governing  system  of  equations  are  rather  general, 
and  are  those  permitted  by  the  numerical  technique  of  Patankar  and  Spalding  (Ref.  1). 
In  addition,  boundary  conditions  on  e  are  required  which  are  taken  as  e  =  0  at  the  wall 
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and  some  prescribed  free-stream  turbulence  level.  If  the  turbulence  is  taken  as  e  =  0  across 
the  boundary  layer,  the  governing  equations  describe  laminar  flow.  If  the  turbulence 
anywhere  across  the  boundary  layer  is  e  >  0,  regardless  of  how  small,  the  flow  is  turbulent. 
That  is,  the  degree  of  turbulence  here  is  determined  by  the  magnitude  of  "e,  and  no 
distinction  is  made  between  fully  developed  turbulent  flow,  transition  flow,  etc. 

It  should  be  pointed  out  that  the  calculation  of  mass  flux  into  the  boundary  layer, 
entrainment,  is  not  handled  as  suggested  in  Ref.  1  for  turbulent  flow,  because  the 
mixing-length  approach  is  not  followed.  Rather,  the  same  subroutine  as  utilized  in  Ref. 
2  for  laminar  flow  was  used  here  for  both  laminar  and  turbulent  flow.  This  permitted 
extremely  small  to  extremely  large  turbulent  intensities  to  be  handled  using  the  same 
entrainment  calculation  scheme. 

2.4  COMPARISONS  WITH  EXPERIMENTAL  DATA 

For  convenience,  the  ratio  of  the  Reynolds  number  based  on  total  conditions,  Rcq, 
and  the  local  unit  Reynolds  number,  Re^,  which  is  more  commonly  used,  is  given  in 
Fig.  1  as  a  function  of  Mach  number.  This  ratio  is  based  on  isentropic  flow  and  assumes 
that  /X  ~  T^  The  assumption  ju  ~  T^  was  used  in  the  present  calculations  with  f  =  2/3. 


Figure  1.  Reo,s/R®<,o,s  3  function  of  M^  for  7  =  7/5. 
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However,  it  is  not  required,  and  any  viscosity  law  can  be  used.  Solutions  presented  in 
Ref.  2  for  thick  laminar  boundary  layers  with  large  transverse  temperature  gradients  indicate 
essentially  the  same  results  were  obtained  whether  /x  ~  T  was  used  or  Sutherland's  law 
was  used.  Therefore,  it  was  assumed  here  that  ju  ~  was  adequate,  particularly  since 
IX  is  significant  only  over  a  small  portion  of  a  turbulent  boundary  layer.  Figure  1  is 
convenient  for  determining  Reo^^  because  usually  the  free-stream  unit  Reynolds  number 
is  known  for  a  given  free-stream  Mach  number. 

The  investigation  by  Winter  and  Gaudet  (Ref.  12)  provides  a  large  number  of  subsonic 
and  supersonic  data  taken  in  the  Royal  Aircraft  Establishment  (R.A.E.)  8-ft  by  8-ft  wind 
tunnel.  These  data  are  particularly  useful  for  evaluating  the  present  work  in  that  the 
development  length  was  established  and  numerous  detailed  measurements  were  made  at 
various  Reynolds  numbers.  Winter  and  Gaudet 's  subsonic  velocity  profile  data  at  = 
0.2  for  16.4  X  10^  <  R®«„x  224  x  10^  are  compared  with  calculated  profiles  in  Fig. 

2.  These  data  were  relied  upon  for  determining  C  in  Eq.  (11)  as  given  by  Eq.  (13).  With 
C  taken  as  a  constant  of  3.5,  Cf  varied  from  about  five  percent  above  the  experimental 


Figure  2.  Theoretical  and  experimental  velocity  distributions 
at  M  =  0.2  and  various  Re  x  ■ 

skin  friction  data  of  Ref.  12  at  Re^^x  “  to  about  ten  percent  below  the  experimental 
data  at  Re_^^x  ==  224  x  10^.  Also,  the  calculated  velocity  distribution  with  C  =  3.5  differed 
somewhat  with  the  experimental  velocity  distribution  for  Re^x  224  x  10^.  It  was  found 
that  if  C  were  changed  to  about  2.7,  then  Cf  agreed  with  the  experimental  skin  friction 
data  at  Re^^x  =  224  x  10^,  and  the  calculated  velocity  distribution  agreed  with  the 
experimental  velocity  distribution.  Of  course,  this  agreement  in  velocities  must  be  the 
case  next  to  the  wall  if  the  skin  frictions  are  to  agree;  however,  not  all  models  attempted 
for  the  various  terms  in  Eq.  (11)  gave  this  quantitative  trend.  It  was  found  that  essentially 
the  same  result  was  obtained  whether  C  was  taken  as  2.7  or  concocted  to  be  2.7  when 
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Re^^x  =  224  X  10^.  Therefore,  C  was  made  to  take  on  the  proper  value  (proper  value 
here  means  giving  good  agreement  with  experimental  data)  at  a  smaller  Re^^x 
=  224  X  10^,  and  Eq.  (13)  was  used  in  the  present  calculations  with  a  lower  bound 
of  2.7.  Actually,  not  a  great  deal  of  effort  was  expended  in  evaluating  or  improving  C, 
and  surely  an  improvement  could  be  made.  However,  the  results  provided  by  using  Eq. 
(13)  were  considered  acceptable,  and  no  further  changes  were  made  for  C. 

Comparisons  between  calculated  and  measured  velocity  distributions  for  the  R.A.E. 
tunnel  (Ref.  12)  are  presented  in  Fig.  3  for  =  2.2  and  26.3  x  10^  <  Re^^x  ^  108 
X  10^.  Comparisons  are  also  made  with  Mach  number  distributions  in  Fig.  4  for  these 
supersonic  flow  data.  Considering  that  C  was  evaluated  for  subsonic  flow,  and  Glushko's 
L  came  from  incompressible  flow,  the  agreement  between  calculated  and  measured  data 
is  considered  good.  The  agreement  improves  as  Re^  x  increases,  with  the  maximum 
discrepancy  at  Re^^x  ”  26.3  x  10^. 


Figure  3.  Theoretical  and  experimental  velocity  distributions 
at  M  =  2.2  and  various  Re  « • 

OO  OOf ^ 


For  comparisons  at  higher  Mach  number  the  experimental  data  of  Gates  (Ref.  13) 
at  M^  =  4  are  considered  in  Figs.  5  and  6.  Two  inflection  points  are  evident  in  each 
of  the  experimental  Mach  number  profiles  in  Fig.  6.  Similar  inflection  points  exist  in 
the  calculated  Mach  number  profiles,  although  the  profiles  at  x  =  4.75  ft  are  in  best 
agreement.  The  solutions  were  obtained  for  a  constant  wall  temperature  of  T^/Tq,^  = 
0.9  as  indicated  by  a  temperature  versus  distance  along-the-plate  plot  as  presented  by 
Gates  (Ref.  13).  However,  the  tabulated  profile  data  in  Table  IV  of  Ref.  13  indicate 
Tw/To,^  =  0.94  for  the  data  at  x  =  2.75  ft  in  Figs.  5  and  6,  which  might  cause  some 
discrepancy.  For  the  same  reasons  as  mentioned  in  discussing  the  R.A.E.  M^  =  2.2  data 
of  Winter  and  Gaudet  (Ref.  12),  the  agreement  in  Figs.  5  and  6  is  better  than  expected. 
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Figure  4.  Theoretical  and  experimental  Mach  number  distributions 
at  M  =  2.2  and  various  Re 

CX>  OOf 


Figure  5.  Theoretical  and  experimental  velocity  distributions 
at  M  =4. 

oo 


16 


AEDC-TR-76-62 


3.0  ANALYTICAL  RESULTS  FOR  TURBULENT  FLOW  NEAR  A  WALL 

Of  interest  here  is  that  portion  of  a  turbulent  boundary  layer  which  extends,  in 
a  normal  direction,  from  the  wall  out  into  the  fully  turbulent  region.  Included  are  the 
so-called  viscous  sublayer  and  buffer  layer  which  are  understood  to  extend  over  values 
of  about  0  to  5  and  5  to  35,  respectively. 

This  portion  of  the  boundary  layer  has  been  investigated  previously  and  summaries 
of  these  investigations,  with  the  known  exception  of  the  work  by  Rannie  (Ref.  14),  are 
given  in  the  books  by  Hinze  (Ref.  15)  and  White  (Ref.  16).  The  approach  of  these 
investigators  included  that  of  assuming  velocity  distributions  in  this  region,  developing 
empirical  expressions  for  the  eddy  viscosity,  or  employing  Prandtl's  mixing-length 
hypothesis  with  various  modifications  near  the  wall.  With  the  exception  of  simply  assuming 
a  functional  form  of  the  velocity  distribution,  none  of  the  previous  investigations  lead 
to  closed-form  expressions  for  the  velocity  except  that  of  Spalding  (Ref.  17);  the  same 
expression  was  obtained  later  by  Kleinstein  (Ref.  18).  Spalding  (Ref.  17)  developed  an 
expression  by  satisfying  the  familiar  logarithmic  velocity  distribution  in  the  fully  turbulent 
portion  of  the  boundary  layer  and  then  matched  this  result  to  an  expression,  obtained 
by  inspection,  which  satisfies  Reichardt's  (Ref.  19)  requirement  that  the  eddy  viscosity 
vary  proportionally  to  y  raised  to  at  least  the  third  power  as  y  becomes  small. 
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Whereas  Spaldings  (Ref.  17)  result  for  the  velocity  distribution  is  accurate  and  is 
in  closed-form,  van  Driest  (Ref.  20)  obtained  one  apparently  as  accurate  by  quadrature. 
Surprisingly,  the  latter  is  perhaps  more  frequently  used.  This  seems  to  be  attributable 
to  the  popularity  of  van  Driest 's  (Ref  20)  physical  development  whereby  he  derived  a 
damping  factor  for  Prandtl's  mixing  length  near  a  wall. 


With  regard  to  these  previous  investigations,  the  present  work  is  of  interest  because 
(1)  the  approach  is  not  founded  on  an  eddy-viscosity  or  mixing-length  concept;  (2)  the 
functional  forms  for  the  velocity  and  Reynolds  stress  distributions  are  not  postulated  but 
result  from  the  analysis;  (3)  simple  closed-form  expressions  are  obtained  for  the  velocity, 
Reynolds  stress,  turbulence  production,  and  direct  mean-flow  energy  dissipation;  and  (4) 
the  results  are  in  good  agreement  with  experimental  data.  The  approach  is  to  model  the 
Reynolds  stress  by  relating  it  to  the  local  turbulent  kinetic  energy  of  the  flow.  The 
momentum  equation,  and  the  turbulent  kinetic  energy  equation  developed  in  the  previous 
section,  are  solved  analytically  in  the  region  near  the  wall.  These  results  are  used  to  obtain 
analytical  expressions  for  the  various  mean-flow  turbulent  quantities. 


The  expression  used  for  the  turbulent  kinetic  energy  is  that  given  by  Eq.  (10),  and 
is,  for  two-dimensional  flow 


de  de 

^  P^Ty  = 


p<\iv>\  de 

dV/dy)  dy 


p<uv>  ^  - 


P<U  V>\ 

dU/dyJ 
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T  = 


p<uv> 


(15) 


(16) 


Eq.  (15)  becomes 


jj,  de  de 

P'^  T.  ^  P^d^ 


1  d  /t  de/dy\  (9U  Ce  t 

o'e  (^y\<5U/(9y  /  ^  dy  d\}/dy 


(17) 


It  is  hypothesized  that  e  is  an  explicit  function  of  U  only.  Then,  Eq.  (17)  can  be  written 


(  TT 


1  d  {  de\  (9U  Ce  T 

dy  \  du)  dy  ~  dU/dy 


(18) 


Using  the  zero-pressure  gradient  momentum  equation 


+  pW 


dy 


dr 

dy 


(19) 
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and  taking  Oe 


1,  Eq.  (18)  simplifies  to 


T 


/9<UV> 


Ce  T 

OU/^y)^ 


0 


(20) 


The  dissipation  length  L  in  the  last  term  in  Eq.  (20)  was  taken  from  Glushko’s  report 
(Ref.  8).  Because  L  was  originally  taken  from  experimental  data  and  behaves  somewhat 
as  a  mixing  length  (see  Ref.  8),  this  length  is  eliminated  by  assuming  it  was  originally 
related  to  the  Reynolds  stress  by  -<uv>  ~  (3U/3y)2  Absorbing  the  constant  of 

proportionality  in  the  constant,  C,  Eq.  (20)  becomes 


T 


p<uv>  — 


Ce  T 
<u  v> 


0 


(21) 


To  obtain  closure,  the  Reynolds  stress  is  modeled  by  relating  it  to  the  turbulent 
kinetic  energy  of  the  flow  according  to  Eq.  (6).  An  exact  value  for  the  constant  in  Eq. 
(6)  is  not  important  here;  rather,  the  significance  of  Eq.  (6)  is  that  the  Reynolds  stress 
is  assumed  proportional  to  the  local  turbulent  kinetic  energy.  Equation  (21)  can  then 
be  written 


r  — „  +  a  Ve  =  (Sr  (22) 

dU^ 

where  a  and  /3  are  constants  and  was  written  for  later  convenience. 

Equation  (22)  is  not  restricted  to  the  region  near  the  wall  nor  to  incompressible 
flow.  However,  rather  simple  and  useful  results  can  be  obtained  by  making  two  assumptions, 
(1)  incompressible  flow  "and  (2)  r  is  independent  of  y  and  is  equal  to  rw  The  second 
assumption  is  common  and  not  without  some  justification.  For  example,  consider  Eq. 
(19).  At  the  wall,  with  no-slip  boundary  conditions,  one  has  that  dr/3y  =  0.  Differentiating 
Eq.  (19)  one  obtains 


The  first  two  terms  on  the  right-hand  side  of  Eq.  (23)  are  zero  at  the  wall  because  of 
the  no-slip  boundary  conditions.  The  third  term  is  zero  by  the  continuity  equation. 
Therefore,  the  first  nonzero  term  for  r  -  Tw  must  be  at  least  of  order  y^ ,  and  has  no 
y  or  y2  term  in  zero-pressure  gradient  two-dimensional  flow  as  sometimes  assumed,  e.g. 
Ref.  21.  Using  these  two  assumptions  of  incompressible  flow  and  r  independent  of  y  (r 
=  Tw),  Eq.  (22)  can  be  written 
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where 


and 


d^E 


du^ 


+  E  =  /3 


(24) 


(25) 


(26) 


The  boundary  conditions  used  for  Eq.  (24)  are  E  =  0  at  u"^  =  0  and  dE/du+  = 
0  at  u"^  =  0.  A  boundary  condition  near  the  outer  edge  of  the  layer  is  not  used  because 
Eq.  (24)  is  valid  only  in  the  vicinity  of  the  wall  because  of  the  assumption  of  r  being 
independent  of  y.  The  boundary  condition  dE/du’*’  =  0  at  u”*"  =0  represents  the  behavior 
of  the  numerical  solutions  of  the  previous  section  near  the  wall,  simplifies  the  final  results, 
and  contributed  to  yielding  a  physically  reasonable  turbulent  kinetic  energy  solution  across 
the  entire  boundary  layer.  Final  justification  for  this  boundary  condition  is  based  on  the 
quahty  of  the  results  obtained.  Using  these  boundary  conditions,  the  solution  to  Eq.  (24) 
is 


E  =  -§  [1  -  cos  («u+)]  (27) 

a 

With  T  =  Eq.  (16)  can  be  written  in  dimensionless  variables  as 

^  =  1  -  aE  (28) 

dy 

Using  Eq.  (27),  Eq.  (28)  can  be  written 

^  =  l  -  ^[1  -  cos  (au+)]  (29) 

dy+ 

Equation  (29)  can  be  integrated  directly  to  give  y'*'  as  a  function  of  u+  if  the  magnitude 
of  the  constant  aj3/a^  is  known,  see  e.g.  Gradshteyn  and  Ryzhik  (Ref.  22).  To  this  end, 
consider  the  fact  that  in  the  fully  turbulent  region  one  has  that  apem  ^  “Tw  (where  e^ 
is  the  maximum  or  peak  in  the  turbulent  kinetic  energy  distribution).  Therefore,  in 
dimensionless  variables 

®  =  1  (30) 


20 


AEDC-TR-76-62 


where  is  a  small  positive  quantity.  That  is  a  small  positive  quantity  is  a  consequence 
of  Eq.  (28)  in  conjunction  with  the  fact  that  E  >  0  and  the  physical  reasoning  that 
du+/dy^  >  0  for  a  smooth,  impermeable  wall.  Then,  using  Eq.  (27)  to  obtain  E^  =  2j3/a2, 
gives 


a 


2 


(31) 


and  Eq,  (29)  becomes 


du+ 

dy+ 


1  + 


+ 


1  - 


cos  (au"^) 


(32) 


Using  the  fact  that  0  <  <1,  and  satisfying  the  wall  boundary  condition  of  u"^  = 

0  at  y+  =  0,  Eq.  (32)  can  be  integrated  (Ref.  22)  to  give 


y 


+ 


2 

ea 


(33) 


Consider  the  limiting  result  of  Eq.  (33)  for  small  u"^.  Neglecting  terms  of  third  order 
in  u+,  Eq.  (33)  gives 


=  u"^  (34) 

which  is,  of  course,  the  classical  sublayer  result.  Away  from  the  wall,  the  velocity 
distribution  depends  on  the  values  of  e  and  a.  The  constant  a  can  be  determined  by 
using  Eq.  (27)  to  force  the  peak  in  the  turbulent  kinetic  energy  distribution  to  correspond 
to  experimental  data  or  to  the  numerical  solutions  of  the  previous  section.  It  can  also 
be  determined  by  taking  e  =  0,  which  corresponds  to  the  assumption  of  apem  =  , 

integrating  Eq.  (32),  and  then  determining  the  value  of  a  which  gives  best  agreement 
with  experimentally  measured  velocity  distributions  near  the  wall.  Fortunately,  the  value 
of  a  obtained  by  any  of  these  methods  is  about  the  same,  and  is,  a  =  0.18.  With  the 
constant  a  determined,  e  can  be  obtained  most  easily  by  fitting  Eq.  (33)  with  experimental 
data  at  the  value  of  u+  where  au'*'  =  tt,  in  which  case  e  =  Tr/ay”**  according  to  Eq.  (33) 
for  au"^  =  TT.  A  reasonable  value  of  e  appears  to  be  e  =  1/8,  Therefore,  Eq.  (33)  becomes 


= 


0.09 


tan"^  -  tan  (0.09  u+) 


(35a) 


Note  that  for  u+  >  7r/0.18,  a  change  of  quadrants  must  be  made,  and 
y"*"  =  ^  +  tan'^  tan  (0.09  for  u"^  >  77-/0.18 


(35b) 
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The  result  corresponding  to  e  =  0  is  more  simple  than  the  expressions  above  and  is  given 
by 

y-*-  =  ^  tan  (0.09  u+)  (36) 

Comparisons  of  the  velocity  distributions  according  to  Eqs.  (35)  and  (36)  are 
compared  with  the  experimental  data  of  Lindgren  (Ref.  23)  in  Fig.  7.  Good  agreement 
exists  between  Eq.  (36)  and  the  experimental  data  for  y*"  less  than  about  100,  and  poor 


Ref. 

Lindgren  (Ref.  23) 
Lindgren (t^ef.  23) 
Lindgren  (Ref.  23) 
Lindgren  (Ref.  23) 
Present 
Present 


Figure  7.  Velocity  distributions  according  to  Eqs.  (35)  and  (36)  and  the 
experimental  data  of  Lindgren  (Ref.  23). 


agreement  exists  for  y*"  >  100.  For  Eq.  (35),  good  agreement  is  obtained  up  to  y^  of 
approximately  200,  and  then  poor  agreement  for  y*"  >  200.  Recall  that  Eq.  (36)  requires 
apej„  =  Tw,  whereas  Eq.  (35)  does  not  have  this  restriction.  However,  both  Eqs.  (35) 
and  (36)  require  that  r  be  independent  of  y,  which  is  not  the  case  in  the  outer  part 
of  the  boundary  layer.  The  discrepancy  between  Eq.  (35)  and  the  experimental  data  in 
Fig.  7  for  y+  >  200  is,  therefore,  attributed  to  the  assumption  of  r  =  ,  and  the  discrepancy 

between  Eq.  (36)  and  the  data  for  y+  >  100  is  attributed  to  the  same  assumption  plus 
the  additional  restriction  of  requiring  apem  =  % .  Equations  (35)  and  (36)  do,  however, 
provide  simple  and  accurate  expressions  for  the  velocity  distribution  throughout  the  viscous 
sublayer  and  buffer  layer  and  into  the  fully  turbulent  portion  of  the  boundary  layer. 
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An  expression  for  the  Reynolds  stress  can  be  obtained  by  using  Eq,  (16)  with  r 
=  Tw  and  either  Eq.  (35)  or  (36)  for  the  velocity  distribution.  The  expression  resulting 
from  using  the  more  simple  velocity  distribution,  Eq.  (36),  appears  to  be  sufficiently 
accurate,  and  is  given  by 


-  <uv> 


sin^  [tan'^  (0.09  y"^)] 


(0.09 

1  +  (0.09y‘^)^ 


(37) 


Equation  (37)  is  compared  in  Fig.  8  with  the  experimental  data  of  Laufer  (Ref.  24)  and 
Kdebanoff  (Ref.  25)  and  that  calculated  by  van  Driest  (Ref.  20)  using  his  damping  factor 
concept.  The  low  Reynolds  number  data  of  Laufer  (Ref.  24)  are  slightly  below  the  other 
data  in  Fig.  8  for  y"^  larger  than  about  30.  However,  the  agreement  in  Fig.  8  is  considered 
good. 


Figure  8.  Reynolds  stress  distributions  near  a  wall. 


Schubauer  (Ref.  26)  considered  the  measurements  of  Laufer  (Ref.  24)  and  Klebanoff 
(Ref.  25)  for  investigating  certain  turbulent  processes  in  pipe  and  boundary-layer  flows. 
Schubauer  (Ref.  26)  points  out  that  the  most  important  outcome  of  the  investigations 
of  Refs.  24  and  25  is  the  revelation  that  the  region  near  the  hypothetical  sublayer  is 
one  of  high  activity  for  turbulence.  Schubauer  (Ref.  26)  presents  a  plot  of  the  turbulence 
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production  and  dissipation  measurements  of  Laufer  (Ref.  24)  and  Klebanoff  (Ref.  25). 
This  plot  is  reproduced  in  Fig.  9  along  with  the  present  predictions  of  the  quantities 
measured  using  Eq.  (36)  for  the  velocity  distribution.  The  present  results  for  the  quantities 
considered  by  Schubauer  (Ref.  26)  are,  for  the  turbulence  production 

^  =  i  sin2  [2  tan-'  (0.09  y^)]  = 

U2  dy+  ^ 

and  the  direct  dissipation  of  mean-flow  energy 


0.09 


a  +  (0.09y+)^, 


(38) 


=  cos^  [tan” ^  (0.09  y"*")]  =/ - ^ ^  (39) 

\1  +  (0.09y+)^/ 


Figure  9.  Distributions  of  turbulence  production  and  direct  dissipation 
of  mean -flow  energy. 


Schubauer  (Ref.  26)  points  out  that  the  maximum  rate  of  production  occurs  at  what 
is  normally  considered  the  edge  of  the  sublayer,  and  that  direct  dissipation  and  turbulence 
production  go  on  at  the  same  rate  at  this  point.  According  to  Eq.  (38),  the  location 
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of  the  maximum  rate  of  production  is  y+  =  11.1  which  is  approximately  the  edge  of 
the  sublayer.  This  is  also  the  location  where  the  direct  dissipation  according  to  Eq.  (39) 
is  equal  to  the  production.  Good  agreement  between  Eqs.  (38)  and  (39)  and  experimental 
measurements  is  demonstrated  in  Fig.  9. 

4.0  ANALYTICAL  RESULTS  FOR  VELOCITY-TEMPERATURE  RELATIONS* 

Relations  between  velocity  and  temperature  in  boundary  layers  were  obtained  in  the 
early  1930's  by  Busemann  (Ref.  27)  and  Crocco  (Ref.  28).  These  results  provide  a 
convenient  means  of  expressing  temperature  as  a  function  of  velocity,  and  countless 
applications  of  these  early  works  have  since  been  made  in  the  study  of  gas  dynamics. 
It  is  interesting  to  trace  the  development  of  the  velocity-temperature  relations,  inasmuch 
as  there  seems  to  be  some  discrepancy  among  the  users  of  these  relations  as  to  the  original 
investigators  and  their  specific  assumptions  and  results. 

The  first  relation  between  velocity  and  temperature  appears  to  be  that  of  Busemann 
(Ref,  27),  who  obtained  the  solution 

CpT  +  U^/2  =  const  =  Cp  T^^^  (40) 

for  laminar  flow  and  Pr  =  1.  For  this  case,  of  course,  one  must  have  T^  =  Tq,^.  Shortly 
after  this  result  by  Busemann,  Crocco  (Ref.  28)  presented  the  solution 

CpT  +  U^/2  =  aU  +  const  (41) 

(where  a  is  a  constant)  for  turbulent  flow  and  Prt  =  1 .  It  is  frequently  misconstrued  that 
Crocco's  (Ref.  28)  result,  Eq.  (41),  was  obtained  for  laminar  flow;  however,  Crocco  (Ref. 
28)  explicitly  states  that  his  consideration  is  that  of  the  equations  for  the  turbulent 
boundary  layer.  Crocco  (Ref.  28)  references  Busemann's  (Ref.  27)  work,  and  later  Crocco 
(Ref.  29)  states  that  Busemann's  (Ref.  27)  result  was  actually  the  particular  solution  portion 
of  Eq.  (41)  corresponding,  to  what  is  presently  called,  an  adiabatic  wall.  Later,  Busemann 
(Ref.  30)  obtained  the  same  general  result  as  Crocco  (Ref.  28),  Eq.  (41),  but  for  laminar 
flow  and  Pr  =  1.  Evaluating  the  constants  of  integration  in  Eq.  (41)  by  using  the  wall 
and  free-stream  conditions,  one  obtains 

T  =  T  +(T  -T)ir-(T  -T)T^  (42) 

W  ^  O,oo  ^  O,oo  oo' 


The  work  in  this  section  was  done  in  collaboration  with  Dr,  M.  D,  High,  ARO,  Inc. 
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for  either  laminar  or  turbulent  flow  with  the  appropriate  Prandtl  number  being  unity. 
If  the  wall  is  adiabatic,  then  +  U_^^/2cp,  and  Eq.  (42)  reduces  to  Eq.  (40). 

E.  R.  van  Driest  (Ref.  31)  later  investigated  the  problem  by  considering  a  variable 
Prandtl  number  across  a  turbulent  boundary  layer.  He  obtained  general  expressions  for 
the  temperature  as  a  function  of  velocity,  Reynolds  analogy  factor,  and  recovery  factor. 
These  general  expressions  require  distributions  of  Prandtl  number  and  shear  stress  across 
the  boundary  layer  which  van  Driest  evaluated  by  following  von  Karman's  (Ref.  32)  idea 
of  separating  the  boundary  layer  into  three  regions,  each  using  different  Prandtl  numbers 
and  shear  stress  distributions.  This  analysis  provided  quite  reasonable  results,  illustrated 
the  effect  of  Prandtl  number,  and  showed  that  the  local  total  temperature  must  exceed 
the  free-stream  total  temperature  near  the  edge  of  the  boundary  layer  for  adiabatic  flat 
plate  flow  with  nonunity  Prandtl  number. 

Walz  (Ref.  33)  and  Michel  (Ref.  34)  have  investigated  the  problem  along  the  lines 
of  van  Driest  (Ref.  31).  By  making  certain  approximations  to  the  general  expressions  of 
van  Driest  (Ref.  31),  for  example,  Michel  (Ref.  34)  points  out  that  if  shear  stress  and 
Prandtl  number  were  taken  as  constant,  one  obtains 

T  =  T„  +  (T,,  -  TJ  IT  -  (T^„  -  TJ  (43) 


where  the  adiabatic  wall  temperature  appears  in  the  definition  of  the  recovery  factor 


r 


T  -  T 

aw  oc 

0  ,OQ  O 


(44) 


It  is  pointed  out  by  Schlichting  (Ref.  35)  and  Michel  (Ref.  34),  for  example,  that  Eq. 
(43)  follows  directly  from  Eq.  (42)  simply  by  introducing  Taw  in  place  of  To,<„,  because 
To^_^  in  Eq.  (42)  is  actually  Taw  for  unity  Prandtl  number.  Equation  (43)  is  frequently 
referred  to  as  the  modified  Crocco  law,  and  it  is  the  same  as  Eq.  (15.19)  in  Schlichting 
(Ref.  35). 


A  parameter  frequently  used  for  presenting  temperature  data  is  T  =  (Tq  -  Tw)/(To 
-  Tw).  From  Eq.  (42)  one  obtains 


T  =  T  (45) 

Note  that  Eq.  (45)  is  for  unity  Prandtl  number  because  of  the  same  restrictions  on  Eq, 
(42).  For  an  adiabatic  wall  Tw  =  Taw  =  and  from  Eq.  (40)  Tq  =  To,^  =  Taw  and 
T  is  undefined.  From  Eq.  (43)  for  an  adiabatic  wall,  i.e.,  Taw  Tw  5  one  obtains  the 
quadratic  form 
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which  does  not  require  unity  Prandtl  number.  However,  Eq.  (46)  does  not  conserve  energy 
across  the  boundary  layer  for  nonunity  recovery  factor.  This  is  easily  seen  by  considering 
the  total  energy  integral  equation  (see  for  example  Shapiro  (Ref.  36)),  which  for  an 
adiabatic  flat  plate  flow  can  only  be  satisfied  if  To  takes  on  values  in  the  boundary  layer 
that  are  both  less  than  and  greater  than  Tq^^,  i.e.,  there  must  be  total-temperature 
"overshoot".  It  is  clear  that  Tq  can  never  be  greater  than  Tq,^  according  to  Eq.  (46) 
assuming  u  <  1. 

Meier  (Ref.  37)  recently  made  an  interesting  investigation  of  temperature  distributions 
in  turbulent  boundary  layers  by  considering  a  mixing-length  hypothesis  to  describe  a 
variable  Prandtl  number  through  the  boundary  layer.  A  somewhat  similar  analysis  for  the 
variable  Prandtl  number  was  carried  out  by  Cebeci  (Ref.  38).  A  variable  Prandtl  number 
approach  certainly  seems  warranted  in  view  of  the  rather  extreme  variations  in  Prandtl 
number  which  have  been  observed  experimentally,  e.g.,  see  Refs.  39  and  40.  However, 
Meier's  (Ref.  37)  approach  to  predicting  velocity-temperature  relations  does  not  conserve 
momentum  or  energy,  does  not  satisfy  the  boundary  condition  at  the  edge  of  the  boundary 
layer,  requires  numerical  solutions,  and  involves  nine  initial  variables  and  parameters.  It 
would  appear  that  results  obtained  using,  for  example,  Refs.  37  and  38,  for  a  variable 
Prandtl  number  description  near  the  wall  in  van  Driest 's  (Ref.  31)  general  expressions 
for  a  variable  Prandtl  number  analysis,  would  be  more  appropriate.  Unfortunately,  such 
solutions  would  likely  need  to  be  carried  out  numerically,  which  would  limit  their  practical 
applications  for  such  purposes  as  data  reduction. 

The  objective  of  the  present  work  was  to  obtain  an  analytical  description  of  the 
temperature  as  a  function  of  velocity  throughout  a  turbulent  boundary  layer  for  constant 
but  nonunity  Prandtl  number.  The  approach  was  to  use  the  equation  resulting  from 
combining  the  boundary-layer  momentum  and  energy  equations  and  model  the  local 
shear  stress  by  expressing  it  as  a  function  of  the  local  turbulent  kinetic  energy  in  the 
boundary  layer.  A  second-order,  nonlinear,  ordinary  differential  equation  results,  for  which 
zeroth-  and  first-order  perturbation  solutions  were  obtained  for  temperature  as  a  function 
of  velocity  in  terms  of  the  assumed  small  parameter  e  =  1  -  Pr^.  Crocco's  (Ref.  28) 
result,  Eq.  (42),  is  recovered  identically  by  the  zeroth-order  solution.  A  Reynolds  analogy 
factor  is  also  calculated  from  the  present  solution. 

4.1  DEVELOPMENT  OF  THE  BASIC  EQUATION 

The  starting  point  of  this  analysis  is  the  equation  obtained  by  making  the  proposition 
that  temperature  is  a  function  of  velocity  only,  T  =  T(U),  and  then  combining  the 
boundary-layer  momentum  and  energy  equations  (see,  for  example,  page  342  of  Schlichting 
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(Ref.  35)).  It  is  shown  in  Schlichting  that  the  proposition  is  correct  if  (1)  the  pressure 
and  wall  temperature  gradients  are  zero  and  the  Prandtl  number  unity  or  (2)  the  wall 
is  adiabatic  and  Prandtl  number  unity  for  nonzero  pressure  gradient.  Because  the  present 
approach  is  to  obtain  solutions  for  Prandtl  numbers  near  unity,  it  is  assumed  that  the 
solutions  are  approximately  correct  for  conditions  (1)  and  (2)  above  for  near  unity  Prandtl 
numbers. 


The  boundary-layer  equations  considered  consisted  of  both  the  laminar  and  turbulent 
contributions  to  viscosity  and  thermal  conductivity,  van  Driest  (Ref.  31)  considered  these 
same  equations,  and  obtained  for  constant  mixed  Prandtl  number,  Pr^ ,  the  expression 


(47) 


Equation  (47)  is  vahd  for  dp/dx  =  0,  or  for  Prm  =  1  and  adiabatic  wall  if  dp/dx  i=^ 
0.  The  shear  stress  in  Eq.  (47)  is  comprised  of  both  a  laminar  and  turbulent  portion. 
Again,  the  Reynolds  stress  is  related  to  the  local  turbulent  kinetic  energy,  e,  in  the  boundary 
layer  by  Eq.  (6).  By  fitting  analytical  approximations  to  numerical  solutions  and 
experimental  data  for  e,  for  values  of  y  larger  than  the  y  corresponding  to  the  location 
in  the  boundary  layer  of  maximum  e  (emax)  as  predicted  by  the  numerical  solutions, 
the  expression 

~ —  =  exp  (48) 

ax 

was  obtained  as  an  approximate  analytical  fit  where  rj  =  y/8  and  c  is  a  constant  equal 
to  4.  An  example  of  the  quality  of  fit  is  given  in  Fig.  10  by  comparison  with  the  present 
numerical  solutions  and  the  experimental  data  of  Klebanoff  (Ref.  25).  Using  Eq.  (48), 
the  expression 

T  =  \J^  p  exp  (-077^'^^)  (49) 

was  obtained  for  an  approximation  to  the  total  (laminar  and  turbulent  portions)  shear 
stress  distributions  where  Ut-^  =  /Pw  • 

Using  Eq.  (49)  for  the  shear  stress,  and  considering  the  terms  containing  r  in  Eq. 
(47),  one  obtains 


^  du  du  P  du  du  Ti 


(50) 
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In  order  to  eliminate  r\  from  Eq.  (50),  it  was  assumed  that  u  =  77^/"^  for  this  term.  This 
assumption  is  justified  only  by  the  final  results  (a  value  of  m  =  7  was  used  for  the  present 
results).  Equation  (50)  can  then  be  used  in  Eq.  (47)  to  obtain 


n— 1 


1  dh  jdh 

h  du  / du 


+  Pr^Cy  -  1)m£  =  0 


(51) 


Figure  10.  Theoretical  and  experimental  turbulent 
kinetic  energy  distributions. 

Introducing  e  =  1  -  Pr^  and  the  terms  as  defined  in  the  nomenclature,  Eq.  (51)  can 
be  written 


h^-  h  .  1)^  .  (1  -  .)Ah=  0  (52) 

which  is  the  basic  equation. 

4.2  CONSTANT  WALL  TEMPERATURE  SOLUTION 

Equation  (52)  is  a  second-order,  nonlinear,  ordinary  differential  equation  for  h(u). 
This  equation  was  not  solved  in  closed  form.  Zeroth-  and  first-order  perturbation  solutions 
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were  obtained  in  terms  of  the  small  parameter  e.  The  boundary  conditions  for  the  constant 
wall  temperature  problem  are 


li(o)  =  h  and  h(l)  =  1 
w 

Assuming  a  solution  of  the  form 

h(u)  =  +  €  hi  (u)  +  •  •  * 


(53) 


(54) 


one  obtains  to  zeroth  order 


=  -A 


and  to  first  order 


du^ 


I  = 


7M~1 


du 


The  boundary  conditions  given  by  Eq.  (53)  become 


(55) 


(56) 


1  andhi(o)  =  0,  hjd)  =  0  (57) 

The  solution  to  Eq,  (55)  for  the  zeroth-order  boundary  conditions  in  Eq.  (57)  is 

h„(ir)  =  hju  -  4^2  (58) 

Equation  (58)  was  used  in  Eq.  (56)  to  solve  for  hi(u).  The  integration  was  tedious  but 
straightforward.  The  first-order  boundary  conditions  in  Eq.  (57)  were  satisfied  by  the  two 
integration  constants  from  Eq.  (56),  and  the  solution  hi(u)  was  included  in  Eq.  (54) 
along  with  Eq.  (58)  to  obtain,  up  to  first  order,  the  solution 


h(5)  =  li„  hju  -  y  u2  +  f 


/3  An  (  _ 

(a  +  l)(a  +  2)  \  ^  / 


+  45(1  -  10  +  (1  -  Df(o)  +  sf(l) 
aia  +1)  ^ 


(59) 


where 


f(5)  =  4^^  -  A)l“  I'C  -  A1  -  (i-  +  A)  In  If  +  Al] 
^  -  K-  A5 


(60) 

(61) 
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and 

A  =  [2Ah„ h,)']’"'  (62) 

4.3  ADIABATIC  WALL  SOLUTION 

The  boundary  conditions  for  the  adiabatic  wall  problem  are 

=  0  and  Rl)  =  1  (63) 

u=0 

Assuming  again  the  form  of  the  solution  given  by  Eq.  (54),  the  same  zeroth-  and  first-order 
equations,  Eqs.  (55)  and  (56),  were  obtained.  The  appropriate  boundary  conditions  become 

i(l)  =  0  (64) 

The  solution  to  Eq.  (55)  for  the  zeroth-order  boundary  conditions  in  Eq.  (64)  is 

EJS)  =  1  +  |(1  _  (65) 

As  before,  hRu)  was  obtained  by  using  the  solution  for  ho(u),  Eq.  (65),  in  Eq.  (56) 
and  satisfying  the  first-order  boundary  conditions  in  Eq.  (64).  The  complete  solution, 
up  to  first  order,  is 

=  1  .  A  (X  _  ^2)  ,  +  i  (1-^)  -  «1)  -  fra]  (66) 


dT  •  ir=o 


dhi 


=  0,K  (1)  =  1  and^  _ 

TT-n  °  du  u=t 


0,  h 


dh 

du 


where 

f®  =  A[(^_  A)ln  14-  -  A|  -  (C  +  A)  In  K  +  A|]  (67) 

4-  =  -  Au  (68) 

and 


A  =  (2AH„ 


(69) 


4.4  REYNOLDS  ANALOGY 

The  derivation  of  the  Reynolds  analogy  factor  s  is  straightforward,  and  for  the 
boundary-layer  equations  considered  is  given  by  (Ref.  31) 
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s 


PV  ('■aw 


'■w) 


du 


u=0 


(70) 


Values  of  s  were  obtained  by  using  Eq.  (59)  for  h(u)  in  Eq.  (70).  The  term  dh/duiu=o 
was  found  not  to  approach  zero  at  the  same  rate  as  (haw  ■  Kw)  for  near-adiabatic  wall 
conditions  for  all  Mach  numbers,  and  a  singularity  occurred  at  hw  =  haw  •  However,  quite 
reasonable  values  of  s  for  hw ,  both  greater  than  and  less  than  haw  »  were  obtained  and 
will  be  compared  with  experimental  data. 


4.5  COMPARISONS  WITH  EXPERIMENTAL  DATA 


To  apply  the  analytical  velocity-temperature  results,  it  is  necessary  to  know  the  Prandtl 
number.  Specific  values  of  the  turbulent  Prandtl  number  have  been  suggested  by  van  Driest 
(Ref.  31)  as  Prt  =  0.86,  Bradshaw  (Ref.  41)  as  Prt  =  0.91,  and  Elser  (Ref.  42)  as  Prt 
=  0.921.  However,  as  mentioned,  experiments  (Refs.  39  and  40)  indicated  that  Prt  varies 
from  above  one  near  the  wall  to  about  0.5  or  0.7  at  large  distances  from  the  wall.  These 
data  have  a  relatively  significant  amount  of  scatter  as  might  be  expected  due  to  the 
experimental  difficulties.  Because  the  variation  in  experimental  recovery  factors  is  less  than 
in  experimental  Prandtl  numbers,  it  was  decided  to  determine  Pr^  by  assuming  a  constant 
r  and  calculating  Prm  using  the  adiabatic  wall  solution.  It  was  assumed  that  the  Pr^ 
determined  by  this  procedure  was  also  applicable  to  nonadiabatic  wall  conditions.  This 
calculation  was  carried  out  by  using  Eq.  (66)  with  u  =  0  and  Eq.  (44)  to  provide  two 
equations  for  the  two  unknowns  Pr^  and  haw  (or  Taw).  The  values  of  Pr^  determined 
in  this  way  are  presented  in  Fig.  11  for  r  =  0.87,  0.88,  0.89,  and  0.90.  Available  values 
of  experimental  recovery  factors  seem  to  indicate  a  value  of  r  =  0.88  as  being  representative, 
and  this  value  was  used  for  the  present  results.  For  r  =  0.88,  0.800  <  Prm  <  0.914 
for  all  Mach  numbers.  For  M^  >  2.5,  the  Prm  from  Fig.  11  for  r  =  0.88  is  essentially 
bounded  by  the  suggested  turbulent  Prandtl  number  values  of  van  Driest  (Ref.  31)  and 
Bradshaw  (Ref.  41). 

Some  of  the  velocity-temperature  results  are  presented  in  u  -  T  coordinates.  It  should 
be  pointed  out  that  T  is  a  sensitive  parameter.  For  example,  for  an  adiabatic  surface 
in  moderate  supersonic  flow  T  =  (To/To,„  -  Taw/To,J/(l  -  Taw/To.J,  and  (1  -Taw /To, J 
is  of  the  order  1/10.  Therefore,  if  a  Tq/Tq,^  overshoot  of  one  percent  occurs,  the  result 
is  a  ten  percent,  or  order  of  magnitude  more  overshoot  in  T.  The  adiabatic  wall  solution, 
Eq.  (66),  is  presented  in  Fig.  12  for  7  =  7/5  and  various  M^  (actually,  rather  than  7 
and  M^,  A  =  (7-l)M^2  the  appropriate  parameter).  The  overshoot  in  T  decreases  as 
M^  increases,  a  result  to  be  checked  subsequently  by  comparison  with  experimental  data. 
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Figure  11.  Mixed  Prandtl  number  as  a  function  of  for  7  =  7/5 


Figure  12.  Velocity -temperature  relations 
for  an  adiabatic  wall  according 
to  Eq.  (66). 
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Because  the  present  result  is  an  approximate  analytical  solution  with  constant  Prandtl 
number,  it  was  of  interest  to  compare  these  analytical  results  with  numerical  solutions 
to  the  full  boundary-layer  equations  including  laminar  and  turbulent  terms  with  variable 
Prandtl  number.  This  seemed  particularly  warranted  in  view  of  the  recent  results  of  Meier 
(Ref.  37).  A  comparison  is  presented  in  Fig.  13  between  Eq.  (66)  and  numerical  results 
provided  by  Adams  (Refs.  43  and  44)  for  y  =  7/5,  =  4,  and  adiabatic  wall,  using 

- Numerical  Results  of  Adams  with 

Cebeci's  (Ref.  38)  Variable  Prandtl  Number* 

- —Numerical  Results  of  Adams  with 


Figure  13.  Velocity-temperature  relations  for 

constant  and  variable  Prandtl  number, 
adiabatic  wall,  =  4,  and  y  =  7/5. 

constant  laminar  and  turbulent  Prandtl  numbers  of  0.88  and  also  the  variable  turbulent 
Prandtl  number  description  of  Cebeci  (Ref.  38)  with  Pr  =  0.71.  (Adams'  calculation  scheme 
was  used  because  of  its  variable  Prandtl  number  capability.  The  numerical  scheme  discussed 
in  Section  2.0  was  run  for  constant  Prandtl  numbers  of  0.88  for  comparison  with  Adams' 
results,  and  the  results  were  essentially  identical.)  According  to  Adams*,  the 
inflections  in  the  numerical  results  for  the  variable  Prandtl  number  for  u  >  0.6  were 
sensitive  to  the  modeling  of  the  variable  Prandtl  number.  Because  of  the  rather  significant 
size  of  the  uncertainty  envelope  near  the  wall  of  the  experimental  data  used  to  qualify 
the  variable  Prandtl  number  analysis  of  Cebeci  (Ref.  38)  (attributable  to  experimental 
difficulties  extremely  close  to  the  wall),  and  because  almost  the  entire  boundary  layer 


*Unpublished  data  furnished  by  J.  C.  Adams,  Jr.,  VKF,  AEDC. 
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corresponds  to  u  >  0.6  for  a  turbulent  boundary  layer,  the  agreement  in  Fig.  13  between 
the  more  exact  numerical  results  and  the  approximate  analytical  result,  Eq.  (66),  is 
considered  good.  Included  in  Fig.  13  are  the  linear  and  quadratic  results  given  by  Eqs. 
(45)  and  (46).  In  addition,  an  analog  computer  solution  was  obtained  for  Eq.  (52)  for 
these  conditions,  and  the  results  cannot  be  discerned  from  the  result  given  by  Eq.  (66) 
in  Fig.  13.* 

Experimental  adiabatic  (or  nearly  so)  wall  data  of  Gates  (Ref.  13)  for  ~  4  and 
Voisinet  (Ref.  37)  and  Voisinet**  for  ~  3  are  compared  with  Eq.  (66)  for  M^=  3.5  in 
Fig.  14.  It  should  be  pointed  out  that  Gates'  temperature  data  were  interpolated  for  u 
less  than  about  0.65,  but  were  measured  for  u  >  0.65.  This  limitation  was  due  to  probe 
size.  Voisinet's  data  were  taken  with  a  fine-wire  temperature  probe  (Ref.  45)  and  do  not 
have  this  restriction.  The  wiggle  or  inflection  point  in  these  data  near  u  of  0.6  is  similar 
to  the  predictions  of  Meier  (Ref.  37)  and  Adams  using  a  variable  Prandtl  number  through 
the  boundary  layer.  According  to  Voisinet**,  this  trend  was  repeatable  in  the  experiments; 
however,  the  absolute  value  of  the  data  in  terms  of  T  might  not 

Sym  ^00  Ref.  Model 

o  4. 11  37  Fiat  Plate 

□  3.96  37  Boundary-Layer  Channel 

^  2.94  VoisineF’^  Boundary-Layer  Channel 

o  2.90  VoisIneF’*  Boundary-Layer  Channel 

1.0 

0.8 

T  0.6 
0.4 

0.2 

0 

0  0.2  0.4  0.6  0.8  1.0 

u 

This  analog  computer  solution  was  provided  by  J.  A.  McClure,  ARO,  Inc. 

**Unpublished  data  furnished  by  R.  L.  P.  Voisenet,  Naval  Ordnance  Laboratory,  June  1974. 


Figure  14.  Theoretical  and 
experimental 
velocity-temper¬ 
ature  relations 
for  of 
approximately 
3  to  4  and 
adiabatic  walls. 
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be  precisely  repeatable.  For  example,  a  one-Kelvin-degree  change  in  the  measured  local 
total  temperature  near  u  =  0.6  would  produce  a  15-percent  change  in  T  for  the  conditions 
of  these  data.  Such  sensitivity  should  be  kept  in  mind  when  considering  data  in  u  -  T 
coordinates. 


Experimental  data  similar  to  those  in  Fig.  14,  except  for  being  approximately 
5  and  6,  are  compared  to  Eq.  (66)  for  =  5.5  in  Fig.  15.  The  quality  of  agreement 

between  experiment  and  theory  is  about  the  same  as  obtained  in  Fig.  14,  which  is 
considered  reasonable. 
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Figure  15.  Theoretical  and  experimen¬ 
tal  velocity-temperature 
relations  for  of 
approximately  5  to  6 
and  adiabatic  wails. 


Mention  should  be  made  of  the  data  of  Sturek  and  Danberg  (Refs.  48  and  49)  which 
were  used  by  Meier  (Ref.  37)  for  the  selection  of  constants  which  appeared  in  his  work. 
Sturek  (Ref.  49)  states  that  these  measurements  were  made  in  essentially  adiabatic  flow 
after  25  to  30  minutes  of  run  time,  and  no  heat-transfer  effects  were  expected  to  appear 
in  the  data.  These  data,  however,  were  shown  by  Sturek  (Ref.  50)  to  exhibit  an  extreme 
"dip"  in  T  for  u  ^  0.6.  Except  for  the  downstream  adverse  pressure  gradient  these 
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measurements  correspond  to  conditions  similar  to  those  of  Voisinet,  although  Voisinet's 
measurements  were  made  with  a  fine-wire  probe  (Ref.  45)  and  Sturek's  (Ref 
49)  measurements  were  made  with  a  wedge-type  probe.  Sturek  and  Danberg  (Ref.  48) 
present  temperature  data  to  within  0.0028  in.  of  the  wall.  However,  the  wedge-type  probe 
used  was  stated  (Ref.  49)  to  be  0.010  to  0.015  in.  thick  at  the  tip  where  the  measurement 
was  made.  Therefore,  some  of  the  data  near  the  wall  were  interpolated  or  extrapolated. 
Moreover,  the  location  of  the  dip  in  T  corresponds  to  the  approximate  location  of  wall 
influence  on  probe  measurements  for  this  size  probe  as  determined  by  Allen  (Refs.  51  and 
52).  Because  Sturek's  data  correspond  to  approximately  the  same  flow  conditions  as 
Voisinet's,  e.g.  a  Mach  number  of  3.5  for  Sturek's  data  as  compared  to  2.9  for  Voisinet's 
data,  only  the  latter  are  considered  here. 

Diabatic,  zero-pressure  gradient  data  of  Danberg  (Ref.  53),  Hopkins  and  Keener  (Ref. 
54),  and  measurements  made  by  Adcock  and  Peterson  and  reported  by  Bertram  and  Neal 
(Ref.  55),  are  compared  with  Eq.  (59)  in  Fig.  16.  The  cold-wall  results,  both  calculated 
and  measured,  tend  to  be  more  linear  in  u  -  T  coordinates  and  have  less  or  no  T  overshoot. 
The  calculated  curves  in  Fig.  16,  for  the  two  conditions  indicated,  essentially  bound  the 
curves  which  would  correspond  to  the  conditions  of  the  other  data  in  Fig.  16. 


Figure  16.  Theoretical  and  experimen¬ 
tal  velocity-temperature 
relations  for  of  6.0  to 
6.5,  nonadiabatic  walls, 
and  zero  pressure  gradient. 


Sym  ^OD 

6.37 
6.34 
6.50 
6.50 
6.00 
6.00 


o 

□ 

o 

a 

o 


^W^O,  CD 

a  515 
a  782 

a  268 

a446 
a  38 
a  49 


Ref. 

53 

53 

54 

54 

55 
55 


Model 

Flat  Plate 
Flat  Plate 
Flat  Plate 
Flat  Plate 
Hollow  Cylinder 
Hollow  Grinder 


37 


AEDC-TR-76-62 


As  alluded  to,  increasing  T  overshoot  for  the  same  u  is  predicted  by  the  present 
results  for  decreasing  and  adiabatic  wall.  This  result  is  qualitatively  supported  by  the 
data  of  Winter  and  Gaudet  (Ref.  12)  as  shown  in  Fig.  17. 


Figure  17.  Theoretical  and  experimental 
velocity-temperature  relations 
for  of  0.8  and  2.2  and 
adiabatic  walls. 


An  interesting  comparison  between  predicted  and  measured  temperature  distributions 
as  a  function  of  Mach  number  is  presented  in  Fig.  18.  These  data  were  used  by  Winter 
and  Gaudet  (Ref.  12)  to  obtain  a  reliable  relation  between  local  temperature  and  Mach 
number  for  2.2  <  M^  <  6.0.  The  present  analytical  result  for  M^  =  2.2  is  shown  in 
Fig.  18.  There  is  an  effect  of  M^  on  the  present  result,  causing  it  to  move  toward  the 
modified  Crocco  result  (the  linear  curve)  as  M_^  increases.  However,  this  effect  is  slight 
over  the  M^  range  of  the  data  in  Fig.  18,  and  the  M^  =  2.2  result  provides  a  lower 
bound  and  yet  is  in  relatively  good  agreement  with  the  experimental  data. 

Some  results  on  Reynolds  analogy  according  to  Eq.  (70)  are  given  in  Fig.  19.  The 
singularity  which  occurs  in  Eq.  (70)  at  T^/Taw  =  1  appears  to  influence  the  solution 
only  in  a  region  near  T^^v/Taw  “  1-  This  observation  is  based  on  the  fact  that  if  the 
region  near  T^/Taw  =  1  is  excluded,  then  a  smooth  curve  can  be  faired  through  this 
region  which  asymptotically  joins  the  results  for  T^/Taw  <  1  and  T^/Taw  >  i-  However, 
the  results  for  small  values  of  T^  /Taw  are  of  primary  interest  here,  because  this  particular 
distribution  does  not  appear  to  have  been  previously  predicted. 
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Figure  18.  Theoretical  and  experimental  Mach  number- 
temperature  relations  for  adiabatic  walls. 


Figure  19.  Reynolds  analogy  as  a  function 
of  wall  temperature  according 
to  Eq.  (70). 


fw^aw 


A  comparison  of  Eq.  (70)  with  experimental  data  is  presented  in  Fig.  20.  The  solid 
curve  in  Fig.  20  is  Eq.  (70),  and  the  dashed  extension  is  the  limiting  result  of  Eq.  (70) 
for  large  T^/Taw  With  the  exceptions  of  Eq.  (70)  and  the  data  of  Wilson  (Ref.  69), 
Fig.  20  was  taken  from  the  report  of  Cary  (Ref.  70).  Note  that  Cary  plots  2ch/cf  rather 
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than  the  inverse  as  plotted  in  Fig.  19.  Of  particular  interest  here  is  the  behavior  of  2ch/cf 
for  small  values  of  h^  /Hq  Cary  states  that  there  is  an  indication  from  the  experimental 
data  that  decreasing  the  ratio  of  h^/Ho,^  decreases  2ch/cf,  an  observation  predicted  by 
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Figure  20.  Theoretical  and  experimental  Reynolds  analogy  data. 


the  present  result.  This  observation  is  further  substantiated  by  the  results  of  Kozlov  (Ref. 
68)  and  Wilson  (Ref.  69),  both  of  which  are  presented  in  Fig.  20.  Kozlov's  (Ref.  68) 
results  are  empirical  expressions  for  the  local  skin-friction  coefficient  and  Stanton  number 
which  were  derived  independently  from  fits  with  experimental  data.  Wilson  (Ref.  69)  states 
that  the  theoretical  prediction  of  decreasing  2ch/cf  with  decreasing  h^/Ho,^  has  not  been 
reported.  For  larger  values  of  hw/Ho,^,  the  present  result  is  seen  in  Fig.  20  to  be  about 
a  compromise  of  the  results  of  previous  investigators. 


5.0  DATA  REDUCTION  OF  PITOT  PRESSURE  MEASUREMENTS 

The  results  of  the  previous  two  sections,  and  the  work  of  Allen  and  Tudor  (Ref. 
71)  in  determining  skin  friction  from  pitot  pressure  measurements,  are  incorporated  in 
this  section  in  a  data  reduction  computer  program.  Input  includes  spatial  measurements 
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of  pitot  pressure  taken  normal  to  a  wall,  and  output  includes  Mach  number  distribution, 
velocity  distribution,  temperature  distribution,  various  boundary-layer  parameters,  and  skin 
friction. 


The  calculation  of  Mach  number  distribution  is,  of  course,  not  new  and  is  included 
only  to  be  complete.  Further  usefulness  of  pitot  pressure  data  follows  if  the  temperature 
distribution  is  measured  or  some  relation  between  velocity  and  temperature  is  used.  It 
is  assumed  here  that  only  pitot  pressure  data  are  available,  and  a  relation  between  velocity 
and  temperature  is  required.  Commonly  used  relations  are  either  the  Crocco  (Eq.  (42)) 
or  modified  Crocco  (Eq.  (43))  result.  Crocco's  relation  is  essentially  the  assumption  of 
a  constant  total  temperature  across  the  boundary  layer  with  the  requirement  that  the 
wall  temperature  be  equal  to  the  free-stream  total  temperature.  The  modified  Crocco  result 
does  not  require  a  constant  total  temperature;  however,  it  does  not  properly  predict  the 
total-temperature  distribution  for  an  adiabatic  wall  (as  discussed  in  Section  4.).  The 
velocity-temperature  relations  used  here  are  those  obtained  in  the  previous  section.  Using 
a  velocity-temperature  relation,  it  is  possible  to  calculate  velocity  and  temperature 
distributions,  boundary-layer  parameters,  and  skin  friction.  A  listing  of  the  computer 
program  is  given  in  Appendix  B. 

5.1  MACH  NUMBER,  VELOCITY,  AND  TEMPERATURE  DISTRIBUTIONS 


The  data  reduction  program  was  written  such  that  if  the  free-stream  flow  is  subsonic, 
then  the  static  pressure  (which  is  usually  measured  to  determine  M^)  is  input  for  the 
y  =  0  point  (at  the  wall)  and  M^  is  calculated  from  the  perfect  gas,  isentropic  relation 
(Ref.  72) 


(71) 


If  the  free-stream  is  supersonic,  then  M^  is  obtained  by  solving  the  adiabatic,  perfect  gas 
equation  (Ref.  72) 


Pq  ,oo 

■  (y+ 

y-1 

y  +  1 

Po,00 

_(y-  +  2_ 

-  (y  -  1)_ 

(72) 


for  M^  by  the  Newton-Raphson  method  (Ref.  73). 


The  Mach  number  distribution  across  the  boundary  layer  is  obtained  by  solving  the 
Rayleigh  pitot  equation  (Ref.  72) 
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^  _  py  +  DM^I  r  y  +  1  (73) 

Poo  L  2  J  \_2yu'^  -  (y  -  1)J 

for  M  for 


and 


—  >  ((y  +  l)/2)>'/(y“» 


(74) 


for  M  for  Po/p^  <  ((7  +  The  static  pressure  is  assumed  constant  across  the 

boundary  layer,  and  is  input  for  <  1  and  calculated  from  Eq.  (71)  for  >  1 .  Equation 
(73)  is  solved  for  M  by  the  Newton-Raphson  method  (Ref.  73). 


The  Mach  number  distribution  provides  one  equation  for  velocity  and  temperature 
at  each  measured  point  (spatial  position)  in  the  boundary  layer  according  to  the  expression 


(75) 


which,  of  course,  is  a  consequence  of  the  definition  of  Mach  number.  The  left-hand  side 
of  Eq.  (75)  is  known  from  the  Mach  number  distribution,  and  therefore,  U/U^  and  T/T_^ 
can  be  determined  if  another  independent  equation  is  available  to  provide  two  equations 
for  the  two  unknowns,  U/U^  and  T/T^.  The  expression  used  for  this  equation  is  that 
given  by  Eq.  (59)  or  (66).  Equation  (75)  and  Eq.  (59)  or  (66),  whichever  is  appropriate, 
are  solved  simultaneously  for  the  velocity  and  temperature. 

Having  calculated  spatial  distributions  of  velocity  and  temperature,  and  using  p/p^ 
=  T_^/T  because  of  the  assumption  that  the  static  pressure  in  the  transverse  direction  is 
constant,  various  boundary-layer  parameters  such  as  displacement  and  momentum  thickness 
can  be  calculated  by  numerical  integration.  Because  of  experimental  difficulties  associated 
with  making  measurements  extremely  close  to  a  wall,  Eq.  (35)  of  Section  3.0  is  used 
to  describe  the  velocity  distribution  (from  which  temperature  is  determined  using  the  results 
of  Section  4.0)  for  <  140. 
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5.2  SKIN  FRICTION 


The  skin  friction  is  determined  by  the  technique  discussed  by  Allen  and  Tudor  (Ref. 
71).  Allen  and  Tudor  present  plots  of  U/U^  versus  (Moo/Mw)  R®oo,y  ^i^h  several  curves, 
each  for  constant  Cf.  Each  plot  is  for  different  M^.  These  plots  are  based  on  the 
compressible  law-of-the-wall  as  presented  by  Fenter  and  Stalmach  (Ref.  74).  This  equation 
is  given  by 


5.75  logjQ 


Re 


oo  ,y 


+  5.1 


(76) 


The  technique  of  Allen  and  Tudor  (Ref.  71),  which  is  an  idea  first  proposed  and 
used  by  Clauser  (Ref.  75)  for  incompressible  flow,  is  to  plot  U/U^  versus  (ju^/Mw )  Re^^y 
on  transparent  paper  and  overlay  it  on  the  plot  in  their  report  for  the  appropriate  M^, 
and  determine  Cf  by  the  curve  in  best  agreement  with  the  experimental  data.  Agreement 
cannot  be  expected  all  the  way  across  the  boundary  layer,  i.e.  for  all  values  of  (m«,/mw) 
Re^^y.  This  is  because  the  law-of-the-wall  expression  is  only  valid  in  a  particular  region 
near  the  wall.  The  reader  is  referred  to  Ref.  76  for  an  example  of  the  agreement  sufficient 
to  determine  Cf  by  this  method.  The  trend  is  for  the  experimental  data  to  follow  a  curve 
which  is  parallel  to  a  constant  Cf  curve.  The  value  of  Cf  which  corresponds  to  the  curve 
fit  that  passes  through  the  experimental  data  in  this  portion  of  the  curve  is  the  appropriate 
value  of  Cf  for  the  experimental  data.  That  is,  there  are  data  points  in  a  profile  set  which 
have  the  same  values  of  Cf.  This  technique  has  the  important  practical  advantage  of  not 
requiring  data  from  the  sublayer. 

In  heu  of  plotting  the  data,  skin  friction  is  determined  in  the  present  data  reduction 
program  by  solving  Eq.  (76)  for  Cf  by  the  Newton-Raphson  method  (Ref.  73)  at  each 
point  across  the  boundary  layer.  Assuming  the  viscosity  is  proportional  to  temperature 
to  some  power,  everything  in  Eq.  (76)  is  known  except  Cf.  The  power  law  ju  ~ 
was  used  in  .the  present  data  reduction  because  it  was  that  used  by  Allen  (Ref.  76).  For 
the  experimental  data  considered  thus  far,  the  constant  values  of  Cf  for  each  set  of  data 
(i.e.  the  data  for  each  profile  set)  have  also  been  equal  to  the  minimum  Cf  in  a  profile 
set.  Equal  here  means  within  the  scatter  of  the  data.  This  facilitates  the  selection  of  the 
appropriate  Cf  in  the  computed  set  of  Cf's. 

It  was  pointed  out  by  Allen  and  Tudor  (Ref.  71)  that  the  law-of-the-wall  given  in 
Ref.  74  has  only  been  verified  for  an  adiabatic  or  near-adiabatic  wall.  Also,  Allen  and 
Tudor  further  state  that  the  surface  should  be  smooth  in  the  aerodynamic  sense,  and 
caution  should  be  employed  in  using  the  technique  in  flows  with  large  pressure  gradients. 
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Verification  of  this  technique  was  investigated  by  considering  the  experimental  data 
of  Winter  and  Gaudet  (Ref.  12).  Winter  and  Gaudet  measured  skin  friction  using  a  balance. 
They  also  made  boundary-layer  pitot  surveys  which  were  used  to  determine  U/U^  which 
in  turn  were  used  to  determine  Cf  by  the  present  data  reduction  program  for  comparison 
with  Cf  as  determined  by  the  balance  measurements.  These  results  are  presented  in  Fig. 
21  and  good  agreement  is  obtained  between  these  two  methods  of  determining  Cf. 


Figure  21.  Skin-friction  coefficients  as  determined  by  balance  measurements 
and  data  reduction  of  boundary-layer  pitot  pressure  measurements. 

6.0  TURBULENT  BOUNDARY-LAYER  MEASUREMENTS 
FROM  TUNNELS  16S,  16T,  AND  ART 

6.1  MACH  NUMBER,  VELOCITY,  AND  TEMPERATURE  DISTRIBUTIONS 

Boundary-layer  measurements  have  been  made  on  the  walls,  ceiling,  and  floor  of  the 
16-ft  Supersonic  Wind  Tunnel  (16S)  in  PWT  by  Baker  and  Pate  (Ref.  77)  and  Maxwell 
and  Hartley  (Ref.  78).  These  measurements  were  made  using  a  pitot  pressure  rake  where 
the  spacing  between  probe  centers  was  0.5  in.  near  the  wall  and  1.0  in.  away  from  the 
wall.  Recently,  measurements  have  been  made  of  the  floor  boundary  layer  in  16S  using 
a  traversing  probe  mechanism  to  make  pitot  pressure  and  total-temperature  measurements 
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simultaneously.  Results  from  the  pitot  pressure  measurements  are  presented  in  Figs.  22 
through  24.  These  results  are  for  Station  -2.9  which  is  61.6  ft  (axial  distance)  downstream 
of  the  geometric  nozzle  throat.  The  numerical  solutions  were  started  at  the  nozzle  throat 
and  the  pressure  gradient  corresponding  to  that  appropriate  for  the  flexible  sidewalls  rather 
than  the  straight  floor  or  ceiling  wall  was  used  because  it  was  known  from  the  nozzle 
design  criteria.  The  initial  conditions  used  to  begin  the  numerical  computations  at  the 
throat  were  estimated  by  obtaining  solutions  (beginning  about  75  ft  upstream  of  the  throat) 
through  the  converging  portion  of  the  nozzle  using  the  pressure  distribution  from  assuming 
one-dimensional  flow  through  the  geometric  area  distribution.  The  numerical  results  for 
Station  -2.9  were  relatively  insensitive  to  the  initial  conditions,  so  long  as  the  initial 
boundary-layer  displacement  thickness  did  not  exceed  that  predicted  by  the  numerical 
solutions  up  to  throat. 
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Figure  22.  Theoretical  and  experimental  velocity  distributions 
^  in  Tunnel  16S  for  =  1.6. 


The  agreement  between  calculations  and  experimental  data  in  Figs.  22  through  24 
is  considered  good.  The  maximum  discrepancy  occurs  for  one  of  the  lowest  Reynolds 
number,  =  1.6  profiles  in  Fig.  22  indicated  by  the  circle  symbols.  This  profile 


45 


AEDC-TR-76-62 


Sym  ^OD  Probe  Size,  in.  Source 

o  2.20  0,739x  10^  a  062  16S  Data 

□  2.20  0.718X#  a  032  16S  Data 

o  2.20  1.210x10^  a  062  16S  Data 

^  2.20  1.679x  10^  a  062  16S  Data 

A  2.20  1.679x  10^  a  032  16S  Data 

- 2.20  OJOQxVp  —  Present  Calculation 


U/Ua, 

Figure  23.  Theoretical  and  experimental  velocity  distributions 
in  Tunnel  16S  for  =  2.2. 
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Figure  24.  Theoretical  and  experimental  Mach  number  distributions 
in  Tunnel  16S  for  =  1.6  and  2.2. 
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measurement  was  repeated  one  week  later  using  a  different  pitot  probe,  and  better 
agreement  was  obtained  as  indicated  by  the  square  symbols  in  Fig.  22.  Although  the  digital 
voltmeter  output  of  the  pitot  pressure  was  not  stable  while  the  circle  symbol  profile  was 
taken,  the  discrepancy  is  not  necessarily  attributed  to  probe  vibration  or  transducer 
problems.  Maxwell  and  Hartley  (Ref.  78)  observed  some  time  ago  that  at  low  Mach  number 
the  boundary  layer  varied  with  time.  This  anomaly  was  not  resolved  in  their  work  and 
it  is  not  resolved  herein. 

Boundary-layer  measurements  have  recently  been  made  on  the  floor  and  east  side 
wall  of  the  16-ft  Transonic  Wind  Tunnel  (16T)  in  PWT  using  this  same  traversing  probe 
mechanism.  These  measurements  were  made  at  Station  -8  on  the  solid  walls  between  the 
first  upstream  hole  patterns  at  the  beginning  of  the  transition  region  which  separates  the 
sohd  and  porous  tunnel  walls.  Measurements  on  the  floor  centerline  were  made  first  and 
compared  with  the  present  numerical  solutions.  The  pressure  gradient  used  for  the 
calculations  was  obtained  from  static  pressure  measurements  extending  32  ft  upstream 
of  the  point  of  measurement,  and  the  pressure  distribution  upstream  of  the  static 
measurements  was  obtained  from  solutions  using  the  Potential  Flow  Computer  Program 
(Ref.  79)  and  were  provided  by  Palko,  Todd,  and  Lutz  of  AEDC.  Examples  of  these 
boundary-layer  results  are  presented  in  Fig.  25.  All  measured  boundary  layers  were  thinner 
than  predicted.  The  results  from  the  potential  flow  solution  for  M^  =  0.6  and  0.8  indicated 
a  pressure  difference  between  floor  (or  top)  and  sidewall  due  to  the  difference  in 
contraction  of  the  top  and  bottom  and  the  sidewalls.  (The  top  and  bottom  wall  contraction 
is  the  same  and  each  sidewall  contraction  is  the  same.)  The  sidewall  pressure  was  larger 
than  the  top  or  bottom  wall  pressure  upstream  of  Station  -60,  and  then  the  bottom  and 
top  wall  pressure  was  larger  than  the  sidewall  pressure  downstream  of  Station  -60,  with 
all  pressures  becoming  equal  at  Station  0.  Considering  the  possibility  that  mass  flow  from 
the  floor  boundary  layer  could  be  taking  place  because  of  this  pressure  difference 
downstream  of  Station  -60,  the  probe  was  moved  to  two  feet  below  the  east  wall  centerline 
at  Station  -8  and  further  measurements  were  made.  Examples  of  these  results  are  given 
in  Figs.  26  and  27.  The  predicted  boundary  layers  are  only  slightly  thicker  than  the 
measured  and  better  agreement  is  obtained  than  with  the  floor  data. 

To  investigate  flow  behavior  further  in  16T,  the  east  half  of  the  floor  and  bottom 
half  of  the  east  wall  were  coated  with  oil  to  visually  study  flow  direction  near  the  wall 
for  a  M^  =  0.7  flow  condition.  Two  primary  results  of  this  investigation  should  be  noted. 
First,  oil  streaks  downstream  of  Station  -60  and  upstream  of  Station  -8  definitely  indicated 
a  flow  direction  across  the  bottom  wall,  away  from  the  centerline,  and  up  the  sidewall. 
Secondly,  although  no  pumping  is  applied  to  the  porous  walls  at  M^  =  0.7,  oil  streaks 
indicated  that  some  of  the  holes  were  actually  pumping  in  the  transition  region,  and  also 
some  streaks  indicated  the  flow  was  influenced  by  this  local  pumping  a  few  feet 
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(approximately  two  to  three)  upstream  of  the  transition  region.  Therefore,  in  addition 
to  mass  being  moved  from  the  floor  to  the  sidewall,  it  could  also  be  removed  from  the 
solid  wall  region  existing  between  the  sawtooth  pattern  of  holes  at  the  beginning  of  the 
porous  wall  transition  region,  thereby  thinning  the  boundary  layer. 
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25.  Theoretical  and  experimental  velocity  distributions  on  the  floor 
centerline  in  Tunnel  16T  for  M  =  0.6  and  0.8. 


Results  from  the  total-temperature  measurements  which  were  made  simultaneously 
with  the  pitot  pressure  in  16S  are  presented  in  Figs.  28  and  29  in  u  -  T  coordinates. 
One  set  of  the  =  2.2  data  is  presented  in  Fig.  30  in  terms  of  T/T^  versus  U/U^. 
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Also  presented  in  Fig.  30  are  the  present  results  according  to  Eq.  (59)  and  Crocco’s  result, 
Eq.  (42),  which  assumes  constant  total  temperature  across  the  boundary  layer  and  is 
frequently  used  to  reduce  pitot  pressure  data.  The  sensitivity  of  the  u  -  T  coordinates 
is  again  made  clear  by  the  fact  that  although  the  trend  of  the  present  analytical  results 
is  in  agreement  with  the  experimental  data  in  Figs.  28  and  29,  the  absolute  value  is  not 
in  agreement  for  all  U/U^  for  these  particular  data.  However,  good  absolute  value  agreement 
is  obtained  in  Fig.  30  in  terms  of  static  temperature.  Fortunately,  for  data  reduction 
purposes,  it  is  the  agreement  in  static  temperature  as  a  function  of  velocity  that  is 
important. 
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Figure  26.  Theoretical  and  experimental  velocity  distributions  2  ft  below  east-wall 
centerline  in  Tunnel  16T  for  =  0.7. 

Comparisons  between  experimental  data  and  the  numerical  calculations  of  Section 
2.0  are  presented  in  Fig.  31  for  =1.6  and  2.2.  The  agreement  is  considered  reasonable. 

The  total-temperature  probe  used  to  obtain  the  temperature  results  in  Figs.  28  through 
31  was  0.046  in.  in  diameter  and  constructed  according  to  the  description  given  in  Ref. 
80.  The  center  of  the  total-temperature  probe  was  positioned  to  be  the  same  distance 
from  the  wall  as  the  center  of  the  pitot  pressure  probe  and  was  traversed  simultaneously 
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with  the  pitot  probe.  The  probe  recovery  factor  was  determined  by  relating  the  measured 
probe  temperature  outside  the  boundary  layer  to  the  tunnel  total  temperature  as  measured 
upstream  in  the  stagnation  chamber.  Tunnel  stratification  was  neglected.  The  wall 
temperature  was  measured  by  determining  the  surface  temperature  of  two  0.25-in.  Gardon 
gages  (Ref.  81)  located  on  the  tunnel  floor  near  the  total-temperature  probe.  These  two 
measured  wall  temperatures  were  within  0.2°F. 
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Figure  27.  Theoretical  and  experimental  velocity  distributions  2  ft 
below  east-wall  centerline  in  Tunnel  16T  for  =  0.9. 
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Figure  28.  Theoretical  and  experimental 
velocity-temperature  relations 
in  Tunnel  16S  for  =  1.6. 
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Figure  29.  Theoretical  and  experimental 
velocity-temperature  relations 
in  Tunnel  16S  for  =  2.2. 
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Figure  31.  Theoretical  and  experimental  spatial  distributions  of  static 
temperature  in  Tunnel  16S  for  =  1.6  and  2.2. 
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6.2  SKIN  FRICTION 

Of  particular  interest  here  are  the  tunnel  wall  skin-friction  coefficients  in  16S.  Tunnel 
wall  skin-friction  coefficients  are  required  in  the  boundary-layer  transition  correlation 
method  of  Pate  and  Schueler  (Ref.  82).  Boundary-layer  pitot  pressure  measurements  were 
made  by  Baker  and  Pate  (Ref,  77)  on  the  ceiling  and  east  wall  of  16S  at  Station  5.4. 
These  measurements  were  used  to  determine  Cf  and  5*  which  were  used  in  the  correlation 
of  the  boundary-layer  transition  measurements  made  on  a  12-in.-diam  hollow  cylinder 
located  in  the  16S  test  section.  Skin  friction  was  determined  for  use  in  Ref.  82  by  the 
expression  Cf  =  20/j2j.  where  was  taken  as  the  axial  distance  from  the  nozzle  throat 
to  the  point  of  measurement.  Baker's  measurements  (Ref.  77)  have  been  reduced  using 
the  present  data  reduction  technique,  and  these  results  are  presented  in  Fig.  32.  Also 
presented  in  Fig.  32  are  measurements  at  Station  -2.9  corresponding  to  the  16S  data 
already  presented. 

The  skin-friction  coefficients  used  by  Pate  and  Schueler  (Ref.  82)  (they  used 
=  3  data  only)  were  those  determined  from  Baker's  ceiling  (straight  wall)  measurements 
using  Cf  =  20/^r.  These  results  are  about  25  percent  larger  than  the  =  3  ceiling  data 
in  Fig.  32  using  the  present  data  reduction  technique.  Also,  there  is  considerably  more 
scatter  than  in  the  data  of  Fig.  32.  Using  Cf  =  20/£j,  the  ceiling  data  (Ref.  77) 
produce  larger  Cf’s  than  the  east-wall  data.  However,  the  east-wall  data  have  larger  skin 
frictions  than  the  ceiling  data  using  the  present  technique  as  shown  in  Fig.  32.  The  skin 
frictions  used  by  Pate  and  Schueler  were  about  10  percent  larger  than  the  east  wall  results 
in  Fig.  32.  The  difference  in  Cf  between  ceiling  and  east  wall  data  for  the  same  Mach 
number  is  attributed  to  the  difference  in  upstream  pressure  gradients. 

Skin-friction  measurements  (i.e.  deduced  from  pitot  pressure  measurements  using  the 
present  data  reduction  technique)  have  also  been  made  in  the  Acoustic  Research  Tunnel 
(ART)  and  Tunnel  16T  in  PWT.  These  data,  along  with  those  from  16S,  are  presented 
in  Fig.  33  using  the  correlation  technique  of  Winter  and  Gaudet  (Ref.  12).  These  data 
cover  the  Mach  number  range  from  0.5  to  3.0.  The  east  wall  (Tunnel  16S)  data  at  Station 
5.4  are  slightly  high  in  Fc  Cf  and/or  Fg  Re^.  This  deviation  is  within  the  scatter  of 
the  data  used  by  Winter  and  Gaudet  (Ref.  1 2),  and  the  correlation  is  considered  reasonable. 
Skin  frictions  obtained  by  using  Cf  =  20 /fir  would  have  larger  values  of  Fc  Cf  than  the 
Tunnel  16S  data  in  Fig.  33,  and  acceptable  correlation  would  not  be  obtained. 

In  practice,  skin  frictions  are  frequently  required  for  axial  stations  or  flow  conditions 
where  data  are  not  available.  To  obtain  Cf  for  different  axial  locations  and  flow  conditions, 
the  following  suggestions  are  given.  A  program  is  available  to  calculate  Cf  according  to 
the  theory  of  White  and  Christoph  (Ref.  83).  It  gives  reasonable  results  (Fig.  32)  and 
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can  be  used  to  predict  the  trend  and  hence  extrapolate  or  interpolate  quite  accurately. 
The  same  can  be  said  of  the  present  finite-difference  calculations,  of  which  some 
comparisons  with  data  are  given  in  Fig.  32.  Sufficiently  accurate  interpolations  or 
extrapolations  might  be  made  using  Cf  proportional  to  some  power  of  Reynolds  number, 
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Figure  32.  Tunnel  wall  skin-friction  coefficients 
in  Tunnel  16S  for  M  of  1.6  to  3.0. 


such  as  the  faired  curves  in  Fig.  32.  Another  possibility  is  to  use  the  correlations  of  Winter 
and  Gaudet  (Ref.  12).  Besides  Fc  Cf  versus  FgRe^,  they  also  present  Fc  Cf  versus 
(F5/Fc)Rex.  These  two  plots  provide  a  direct  relation  between  Re^  and  Rex  and  hence 
interpolations  or  extrapolations  can  be  made  graphically. 
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Figure  33.  Correlation  of  tunnel  wall  skin-friction  data  from  three  tunnels  in  PWT. 


6.3  BOUNDARY-LAYER  PARAMETERS 

The  frequently  used  boundary -layer  parameters  of  displacement  thickness,  6*, 
momentum  thickness,  6,  and  shape  factor,  H  =  6*/0,  are  presented  in  Fig.  34  for  Tunnel 
16T  at  Station  -8.  These  data  indicate,  as  did  the  velocity  data,  that  the  floor  boundary 
layer  is  thinner  than  the  east  wall  boundary  layer.  For  example,  the  =  0.6  and  0.8 
floor  data  have  a  5*  of  about  0.47  in.  at  Re^/ft  of  approximately  4.5  x  10^,  whereas 
the  =  0.7  east  wall  data  show  6*  to  be  0.57  in.,  or  about  20  percent  larger.  Although 
the  calculated  shape  factors  are  in  agreement  with  measured  shape  factors  for  both  floor 
and  east  wall  data,  only  relatively  good  agreement  with  5*  and  6  is  obtained  with  the 
east  wall  data.  Plausible  explanations  for  the  difference  in  floor  and  east  wall  boundary 
layers  have  been  given  when  the  velocity  distributions  were  considered. 

Boundary-layer  parameters  in  Tunnel  16S  at  Station  -2.9  are  considered  in  Fig.  35. 
Good  agreement  between  calculated  and  measured  data  is  obtained.  Results  for  one  of 
the  low  Reynolds  number,  =1.6  condition,  are  shown  to  be  low  as  compared  to 
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the  other  data  in  Fig.  35.  This  anomaly  was  also  discussed  earlier  when  the  velocity  data 
were  considered.  Repeat  measurements  of  this  flow  condition  provided  results  in  better 
agreement  with  the  calculations  and  the  other  data  as  shown  in  Fig.  35.  Boundary-layer 
parameters  from  the  measurements  of  Baker  can  be  found  in  Ref.  77.  Boundary-layer 
parameters  for  a  wider  range  of  Mach  numbers  have  been  measured  and  reported  by 
Maxwell  and  Hartley  (Ref.  78). 


Figure  34.  Boundary-layer  parameters  in  Tunnel  16T  for  Station  -8. 
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lO^  10®  10^ 

Figure  35.  Boundary-layer  parameters  in  Tunnel  16S  for  Station  -2.9. 


7.0  CONCLUSIONS  AND  COMMENTS 

The  numerical  boundary-layer  computations  described  in  Section  2.0  have  provided 
a  useful  tool  for  estimates  of  boundary-layer  problems  associated  with  various  research, 
test,  and  test  facility  development  projects  conducted  in  PWT.  Of  primary  importance 
in  this  work  was  the  assessment  of  the  modeling  for  the  Reynolds  stress  and  the  calculation 
of  the  turbulent  kinetic  energy.  Although  the  particular  modeling  of  the  Reynolds  stress 
used  here  has  been  used  by  previous  investigators,  the  particular  modeling  of  the  turbulent 
kinetic  energy  equation  (on  which  the  Reynolds  stress  depended)  has  not  been  used.  The 
approach  followed  permitted  the  use  of  the  natural  turbulent  kinetic  energy  boundary 
conditions  and  the  computation  of  e  throughout  the  boundary  layer  without  having  to 
handle  such  regions  as  the  laminar  sublayer,  buffer  layer,  logarithmic,  and  wake  regions 
separately. 
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The  analytical  investigation  of  turbulence  near  a  wall  (Section  3.0)  provides  convenient 
expressions  of  describing  the  flow  in  this  region.  Of  interest  is  the  fact  that  these  results 
were  obtained  without  the  use  of  mixing-length  theory,  in  which  the  approach  appears 
to  be  unique.  This  was  made  possible  by  the  form  of  the  turbulent  kinetic  energy  equation 
developed  in  Section  2.0. 

The  analytical  results  obtained  in  Section  4.0,  for  temperature  as  a  function  of  velocity 
throughout  the  boundary  layer,  were  shown  to  be  improvements  over  previously  used 
results.  The  application  of  these  results  to  the  calculation  of  the  quantity  2ch/cf,  indicated 
that  2ch/cf  decreased  for  small  decreasing  values  of  h^/Ho,^.  This  prediction  was 
substantiated  by  comparisons  with  experimental  data. 

The  analytical  results  of  Sections  3.0  and  4.0  were  used  in  Section  5.0  in  a  computer 
program  for  the  data  reduction  of  pitot  pressure  measurements  made  in  turbulent  boundary 
layers.  In  addition  to  Mach  number,  velocity,  temperature,  and  various  boundary-layer 
parameters,  skin  friction  is  determined  from  the  pitot  pressure  measurements.  The 
technique  used  to  determine  Cf  has  the  important  practical  advantage  of  not  requiring 
measurements  near  the  wall  in  the  sub-  or  buffer  layers  and  hence  eases  experimental 
requirements. 

The  recent  pitot  pressure  and  total-temperature  measurements  made  in  Tunnel  16S 
were  in  good  agreement  with  the  present  theoretical  results.  However,  this  was  not  the 
case  with  regard  to  the  solid  wall  pitot  pressure  measurements  made  in  Tunnel  16T.  The 
calculated  boundary  layer  was  consistently  thicker  than  the  measured  boundary  layer, 
particularly  the  floor  boundary  layer.  Potential  flow  solutions  indicate,  however,  that  the 
flow  near  the  tunnel  wall  experiences  a  pressure  drop  from  the  floor  toward  the  sidewall, 
and  oil  streak  experiments  conducted  in  the  lower  east  quadrant  for  60  ft  upstream  of 
the  test  section  at  M  =  0.7  indicated  that  a  crossflow  did  occur  from  the  floor  toward 

oo 

the  east  wall.  The  oil  streak  experiments  further  indicated  that  crossflow  was  induced 
over  the  solid  wall  between  some  of  the  upstream  holes  in  the  transition  section,  which 
extends  10  ft  upstream  of  the  test  section,  thereby  effecting  the  measured  boundary  layer. 

The  skin-friction  coefficients  obtained  by  using  the  present  data  reduction  program 
to  reduce  the  measurements  of  Baker  and  Pate  (Ref.  77)  were  not  in  agreement  with 
those  obtained  by  Pate  and  Schueler  (Ref.  82)  and  used  to  correlate  the  Tunnel  16S 
boundary-layer  transition  data.  However,  preliminary  results  indicate  that  the  effect  of 
using  the  present  skin-friction  coefficients  in  the  correlation  of  Pate  and  Schueler  (Ref. 
82)  is  not  necessarily  to  invalidate  the  correlation  of  the  Tunnel  16S  data  as  it  is  to 
extend  the  correlation  below  the  original  lower  Mach  number  bound  of  three  to  include 
the  remainder  of  the  Tunnel  16S  transition  data  (Ref.  77)  (which  include  Mach  numbers 
down  to  two)  not  correlated  by  Pate  and  Schueler  (Ref.  82). 
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An  effort  was  made  to  make  the  present  computer  results  readily  available  for 
engineering  applications.  The  boundary-layer  computer  program  given  in  Appendix  A  has 
simple  input,  and  unless  desired,  no  dimensional  variables  or  properties  are  input  and  no 
units  are  involved.  The  pitot  pressure  data  reduction  computer  program  given  in  Appendix 
B  is  essentially  self-contained  and  has  been  used  for  online  and  offline  data  reduction.. 
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APPENDIX  A 

BOUNDARY-LAYER  COMPUTER  CODE 

Improvements  to  this  code  that  have  to  do  primarily  with  the  region  near  the  wall  were 
made  subsequent  to  its  submission  for  publication.  These  improvements  hardly  produce 
discernible  changes  from  the  results  given  herein.  The  code  used  for  the  computations  in  this 
report  is  hsted  here  and  the  new  version  is  available  on  request. 
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C  MAIN  A  a 

IMPLICIT  REALMS (A-HtO-2) _ _ _ _ 

COMMON  /GEN/  PEI « AMI  9 ANE 9 DPDX9PREF <3) *PR ( 2) «P (2) «DEN« AMU»XU« XDt XP*  A  3 

2IMEN9IHEAT9Z9TO9INTG/I/N9NPI9NP29NP39NEQ9NPH9KEX9KIN9KASE9KRAD/B/B  A  5 
3ETAtGAMA(2)  #TAUIiTAUEiAJlia)_9JLj£lElj_lNDlt2JiINDEt2)/V/Ut20Q)  iFt292  A  6 
400) 9R(200) 9RHO(200) 90M(200) 9 Y (200 ) /C/SC (200 ) 9 AU (200 ) 9BU (200) 9CU (20  A  ? 

50) 9  A (29200)  98(29200) iC (2  9  200) /D/YR (200) 9UR (2001 9RR (200 ) 9HRI2OO) tXM  A  S 
6  (200)  9PITOT(200) 9  TEMP (200 ) /E/DSTAR (300 )  9XRS(300)  9RWRS(300)  9COSAL(3  A  f 

_ JAQ)  -  -  _ ^ _ 

1/F/ITUR89IPRINT9TURBIN9TURBFS9IDELY 

COMMON  /U  AK9ALM6  _  .  _  A  11 

CALUERRS£T(207  9  2569-»l9l) 

CALLERRSET (20892569-191)  _ 

CALLERRSET (20992569-191) 

_  CALLERRSET  125142564- 111)  _ _ _ 

CALLERRSET (25392569-191) 

CALLERRSET (261 9 2569 -111) 


1  CONTINUE  A  12 

AKsUO^O  A  13 

INTG=0  A  14 

_ CALL-  BEGIN _ _ _ _ _ _A._1S 

AMIaO.D^O  A  18 

_  AME®0»D«0  A  If 

OX®XSTEP»YSTART 

. .  XOuXU^DX  . . - 


GO  TO  3  A  20 

CALL  READY _ _ _ _ A _ 2L- 

TME  FOLLOWING  STATEMENT  IS  TO  KILL  SHOT  IF  PROBLEMS  (LIKE  MERGING)  A  22 

OCCUR  IN  SUBROUTINE  READY  A  23 

IF  (XD96T9XL)  60  TO  13  A  24 

CONTINUE  .  _ _A  25 

00  5  MBI9NP3 

■JXXiKNPa-MLAJ  (NP3 » 0 )  GO-  TP— 4 _ 

CONTINUE 

L®NP3-M  _  _ _ _ _ _ 

OELTAsY(L) « (0.990«0«U(L)/U(NP3) )»(Y(L^l)»y(L) )/((U(L«l)“U(L) )/U(NP 

OSTAR(300)®OELTA 


INTGalNTGtl - A — 24  ^ - 

XD®XU«DX 


_ -CALL  PRE  IXOiOPilXL  . . . . —  -A  J9 

CALL  ENTRN  A  40 

_ CALU  PRE  (XUfOPOX)  _ _ _ _  - .  A  41 

C  CHOICE  OF  FORWARD  STEP  A  27 

_ DX«0»i)5D»0»P£ I /  (RiPlAMl-R  tNP-3).tAMEl _ 

OX®DABS(DX) 

_ IF  (DX*GT®XSTEP«»Y  (NP3)  )QX«XSTEP»Y  tNP3i _ _ _  _ _ 

XOaXU^DX  A  29 

_  _ IF  (KASE«Ea®2)  GO  TO  6  -  _  -  A  42 

IF  (KIN»EQ,1)  CALL  MASS  (XU9XO9AMI)  A  43 

_ IF  (KEX^EQal)  CALL  MASS  LXUiXX)! AME) _ i _ 44 - _ 

CALL  WALL  A  45 

-4- _ CALL  OUTPUT  --  _  _ _ _ - _ A  46 

CALL  PRE  (XO9OPOX)  A  47 

-  -__CALL  CQEEE _ _ _ _ _  _ A_4S_ 
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SETTING  UP  VELOCITIES  AT  A  FREE  BOUNDARY  ~  A  49 

MODIFIED  FOLLOWING  STATEMENT  FOR  INTERNAL  CORE__ELQM _ _ _ _ A_iO_ 

IF  (KEX»EQ«a)  U(NP3)aOSQRT(2tD«0)»OSQRTn»O40®PPO«*( (OAM«l.O*0)/6A  A  51 

«M))  _ _ _ _ : _ 

IF  (KIN®EQ,2)  U(n3OSQRT(U<n»Un)-2«O«0»(XO«XU)<»OPDX/RHO(l) )  A  52 

CALL  SOLVE  ( AUf  8U»CU»U*NP3)  _ _ _ i_„  5J ... 

C  SETTING  UP  VELOCITIES  AT  A  SYMMETRY  LINE  A  54 

IF  {KINtNE«3)  GO  TO  7  _ A  5Sl 

U(1)®U(2)  A  56 

IF  (KRADaEQafl)  UIDa. 750*0 « Ui £)  ♦  _ _ i _ 5I_ 

7  IF  (KEX^EQ.S)  U (NP3) ® ^TSO^O^U INP2) * •250*0»U (NPU  A  58 

IF  (NEQoEQ«l)  GO  TO  14  _ A _59  . 

DO  13  gtslsNPH  A  60 

DO  8  l®2fNP2  . . . IL_j6I._ 

AU(I)«A(U9l)  A  62 

BU{I)®BIJiU_^  _  _ _ A _ 6a_ 

S  CU(I)®C(J9l)  A  64 

DO  9  ISUNP3  _ A _ 6l5_ 

9  SCa)®FU»I)  A  66 

CALL  SOLVE  ( AU^BUjCUfSCiNPS)  _ ^ _ _ A _ ^61... 

DO  10  I®I«NP3  A  68 

10  F(JfI)sSC<l>  _ _ A _ 6JL- 

IF  (KASE®EQ*2)  60  TO  U  A  70 

C  SETTING  UP  WALL  VALUES  OF  F  _ _ A_JLl_ 

IF  (KIN.EQ.UANO«|NDI(J)  •£0.2)  F  (Of  1)  «  M I  •D«0*3ETA*GAMA  (U) )  «F  ( J92)  A  72 

1-  ( 1  sO*0*SETA-GAMA  (U) )  »F  (U®3)  J  «a5D*0/6AMAXaj _ A_  .l^-- 

IF  (KEXtEQtUANO^INDElJ)  »EQ«21  F  ( JsNPS)  «  (  ( 1  •0*0*BETA*6AMA  (U) )  •£  ( J»  A  74 

_  1 NF2 >  »  ( 1  i.D»0*BETA”6AMiLLaU^£JUijti!iPll^)  (iiJ _ A _ IS_ 

C  SETTING  UP  SYMMETRY»LINE  VALUES  OF  F  A  76 

11  IF  (KXN«NE«3)  GO  TO  12  _ _ — _ A _ IZ_ 

F(J«i)®F (Jf2)  A  78 

IF  (KRAD»EQ®0)  F  U»  1)  ®a750*0®F  (U»2)  f»aSO*0«»Fia4JJ _ _ _ A _ 79 

,12  IF  (K£X«EQ»3)  F  (UsNPS)  750*0®F  IU»NP2)  * «250*0»F  (U»NP1 )  A  80 

..13 _ CONTINUE  _ _ _ _ _ _ _ A _ ai__ 

14  XP®XU  A  82 

XU®XD  _ _ A  — 8A  - 

PEI®PEI*DX®  (R<U®AMI«R(NP3)®AME)  A  87 

C  THE  TERMINATION  CONDITION  _  _ _ A _ SJL.. 

IF  (XU»LT<iKL)  60  TO  2  A  89 

_ LgHlNT®lNTG«l _ _ _ 

CALL  READY 

_ CALL  WALL  .....  .  .  _ _ 

CALL  OUTPUT 

STOP  _ : _ A _ 91_ 

END  A  92*» 
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begin 


SUBROUTINE  BEGIN  C  1 

IHPLICIT  REAL®8(A»HfO-Z) 

COMMON  /GEN/  PEI  9  AMI « AMEtDPDXf PREF (2) «PR (2) #P (g) sDENt AMU9XU«XD»XP9  C  2 

iXLtDXeXSTEPtCSALFAtALPHAtXRfREORStGAMsZETAtPPOtTWTOsYSTARTiUSUPiIO  A  4 
2IMEN9lHEATfZ®TOf INTG/I/N9NPlfNP2  NP3#NEQ9NPH9KEX»KIN9KAS£tKRAO/a/B  C  4 
3ETA®G&MAJ2)  9TAUI  #TAU£#  AJI  (2)  #  AJE  (2)  t  INDUg)  ®  INDE  ( 2) /V/U  ( 200)  #f  (2«2  C  5 

400) fR(200) 9RHO(200) 90M(200) ®Y(20Q)  C  S 

l/F/ITURitIPRiNTiTURBINsTURBFSf IDEUY 

C  PROBLEM  SPECIFICATION  C  7 

READ{5®28)KRA0iI0IMENfNEQ9KEX9KINiIHEAT«ITURBaID£LY«-lPaiMI»N  _ 

READ  (5f29)  REORS ® ZET A  9  PR ( 1 ) #6AM » ALPHA® XR® XL  9 USUP® YST ART  * TWTO® XSTE  C  9 

lPfPR«2) iTURBINiTURBFS 

lOIMEN^O  FOR  PLANAR  FLOW  AND  IDIMENsl  FOR  AXISYMMETRIC  FLOW  C  U 

INITIAL  EDGE  OF  BOUNDARY  LAYER  IS  YSTART  C  12 

PREF(1)3PR(1)  C  13 

PREF(2)sPR(2) 

C  APPROXIMATE  CALCULATION  OF  UEOGE  FROM  ONE  DIMENSION  FLOW  RELATIONS  C  14 

XUs0«D*0  C  20 

CALL  PRE  (XUsOPOX) 

UEDGEaOSQRT (iaD^O^ ( 1 9D«0“PPO»® ( (6AM«1 @0^0 ) /GAM) ) ) 

KASE^g  C  18 

IF  (K1N»EQ«1®0R«KEX»EQ®1)  KASEiil  C  19 

NPH»NEO-l  C  21 

NPlsN^l  C  22 

NP2sN«2  C  23 

NP3®N^3  C  24 

C  INITIAL  VELOCITY  PROFILE  C  2S 

Y(1)®O«0D^O  _  C  26 

Un)®O,0D^O  C  27 

DELU^UEUGE/OPLOAT (NP2) 

DO  1  I®29NP3  C  30 

U(l)aU(l-l) «DELU 
Y (l)sYiTART® (U(I)/UEOGE)®«IDELY 

1  CONTINUE  C  34 

IF(ITUR8,EQ»l)60  TO  30 
DO  30  1®2®NP3 
ETA^Y(I)/YSTART 

Ua)«(2#O^0»£TA-ETA»ETA)®UED6£ 

30  CONTINUE 


C  CALCULATION  OF  SLIP  VELOCITIES- AND- OISXAMCXS _ X _ 35— 

BETA^sKO^O 

60  TO  (29394)  9  KIN  C  37 

2  U(2)sU<3)/(l90'i-0*290«0®B£TA)  C  30 

Y(i)®Y(3)®8ETA/(2®O4O«0ETA)  C  39 

60  TO  6  C  40 

3  UilsU(l)®U(U  _ - _ ^ _ C _ 40 _ 

U13^U{l)»U(3)  C  42 

U33®U(3)»U<3)  C  43 

SQ804®D«O»Un«»l2®O«O®U13^9«D^O^U33  C  44 

U(2)s(l6,090®Ull--4«0«0»U13«U33)/<2®0#04«U(l)^U(301^DSQRTISQ)l  C  45 

Y(2)8Y(3)»(U(2)^U(3>»2®D*0»U(1)  )®950«0/(U(2)<-U(3)«U(n  )  C  46 

60  TO  6  -C _ 4I_ 

4  IF  (KRAD«NE,0)  GO  TO  i  C  48 

U(2)«(4«D40«U(n«U(3)  )/3.D40  ,  - _ _ _ _ _ C  49 

Y(2)809D«0  C  50 

60  TO  6  _ 
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5  U(2)®UCl)  C  52 

Y(a)*y<3)/3.D#0  C  53 

60  TO  (79S#9)f  KEX  C  54 

U(NPg)^UCNPl)/(i.U^0^2*D«0»BETA)  C  55 

y<NP2>«Y(NP3)»(Y{NP3)-Y(NPl)  )<»BETA/(2«D«0*BETA>  C  56 

60  TO  10  C  57 

UnsU(NPl)«U(NPU  C  5B 

Ui38U(NPl)«U(NP3)  C  59 

U33sU<NP3)®U(NP3)  C  60 

SQ»s84«.D^Q«U33®i2«0*O®Ul3^9.D*a»Uli  __  - . . . . .  .  C  61 

U(NP2)»U6»D^O«»U33»4»0^0»U13«UU)/(2«0«0»<U(NPU»U(NP3)  )^OSQRT(SQ)  C  62 

Y(NP2)sY<NP3)®(Y(NP3)»Y(NPU  >«(U(NP2) *0 (NPl) »2*O*0*U (NP3) ) •.5D*0/ (  C  63 

lUiNPZ)  ^U(NPn  ♦UCNP3)  )  C  64 

60  TO  10  C  65 

g(NP21  s(4»O«0»U(NP3)»U<NP1)  _ _  .  C  66 

Y(NP2) aY(NP3)  C  67 

CONTINUE  C  68 

IF  (NEQ.EOa)  60  TO  20  C  69 

DO  19  J®itNPH  C  70 

INITIAL  PROFILES  OF  OTHER  DEPENDENT  VARIABLES  C  71 


C0NST»U£06E«UED6£«TURaiN  . . 

C0NSFSsUEDGE#UE06E»TURBFS 
T£DGE®l*O^O^U£D6E»«2/2«0^0 
TAWTO^O  s.8i0#0«  <  KD«0-T£D6£)  ♦TEOGE 
IF{IHEAT®EQ«0)TWTOiiTA«TO 

DO  U  lalfNPS  C  72 

£TA«YCI)/YSTART  . . . . . 

IF<j9£Q®l)F<U,I)sTWTO^(TAWTO»TWTO)«(U(I)/UED6£)^(UD^O-TAWTO)<» 

I  (U(I)/U£06E)»»2 

IF (J»EQ929AND,ETA«LT«0«03D^0)F(J,n®CONST«DERF(72,3O*0*ETA) 
IF(J»EQ«2»ANO«ETA.6E«0,03D«0)F(J«1)bCONST«(1«0«0«OCOS(3«I4|592650« 

10»ETA)  J/2®D«-04CONSFS»(ETA“0.030*01 
_ llLUU£fl«2^NOU^^llF(U«l)»0*D«a  _ _ _ 


CONTINUE  C  75 

CALCULATION  OF  CORRESPONDING  SLIP  VALUES  C  76 

6AHA(J)®1«0*0 

GO  TO  a2iUtl4)  e  KIN  C  78 

12  FUf2)®F(Jf  U^IF(J,3)«F<J9n  l»a®0*O^BETA-GAMA(J)  )/(I.040«BETA«6AM  C  79 

- - »A^UU -  -  -  - • -  - -  -  . 

GO  TO  15  C  80 

13  G8(U(2)^U(3)«S«D<>0«U(l))/(5»0«0«(U(2)«U(3)l-«8.0#a«Uan  C  81 

GF®(UD«0^PR£F(J)  )/(UD«0*PREF(jn  C  82 

GFa(a#6F)/(UD^0^Q<»6F)_  C  83 

F(J®2)*F<J@3)^6F^(1,0*0-GF)»F(J,U  C  84 

_ _ ao-io-  X5- _ c___a5 

14  F(J92)®F<UfU  C  86 

if  CKRADeEQ.O)  F  (U#2)  8(4afi#-04FUf  l)®F4J»a3)/a*O«0  --  C  87 

15  GO  TO  (16«179ie}f  KEX  C  88 

16  F(JfNP21«FU»NP3)^(F(U«NPXl'^F(U»NP3n#(Ua«O^BETA-6AMAUXiy(l^O^O^  C  89 

18£TA^6ANA(U> >  C  90 

_ GO  ja  is _  ■ _  _ _ C  94^_ 

17  6s<U<NP2)«U(NPl>«0,O«O®U(NP3)l/(5.O«O»(U(NP2)*U<NPl))^8.O«O«»UlNP3)  C  92 

GF«a®D^0-PREF(jn/(l«D#O*PREF(J)  )  C  93 

_  Gr»CG#GF)/(l®D^0«6*6F)  _  _  _  C  94 
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BE3IN 


10 

19 

20 

C 


21 

C 

22 

23 

24 
C 


25 


26 


27 

C 

20 

29 


F  <JfNP2)sF<JtNPU*©F^(UD#0-6F)^F(JfNP3) 

60  TO  19 

F (JfNP2)s(4,O^0®F(JfNP3)-F(J«NPl) )/3.D«0 

CONTINUE  .  -  .  _ 

CONTINUE 

CALL  OENSTY  _ _ _  _ 

CALCULATION  OF  RADII 

CALL  RAO  (XUfRa)  fCSALFA)  .  -  - 

IF  «CSALFA«EQ0O®O«OtOR«KRADpEQ*O)  60  TO  22 

DO  21  I^2fNP3  ^  ^  _ _ _ 

R(I)®R(1)-V(I)»CSALFA 

CHANGE  made  IN  STATEMENT  NUMBER  28  FOR  INTERNAL  FLON _  _ 

GO  TO  24 
00  23  I®2@NP3 
Ra)®R(U 

CONTINUE  - _ _ 

CALCULATION  OF  OMEGA  VALUES 

0M(1)  siOeO^O  _  .  -  _ 

OM(2)®O0O^O 
DO  25  IS39NP2 

OM(I)sOMCI®l)  ^95O^0«^(RHO(I)^U<n^R(I)  ^RHO(l-l)®U(I«»n«R(I-l)  )«^(Y(I 

U-Y(l-l))  _ 

PEIsOM(NP2) 

00  26  I®3#NP1  . .  . . . . 

0M(I)*0M(I)/PEI 

0M{NP2)  SU04.0  _ _  -- 

OMINPSlislaO^O 

IF  <NEQ®Ea®l)  RETURN  _ _ 

DO  27  J*lfNPH 

IF  iKEKtEQeU  INDEU)®!  _ 

IF  CKIN^EQaU  INDI(J)®1 

CONTINUE  -  .  _ _ ■  _ 

RETURN 


F0RMAT(@n«2I3) 
FORMAT  (SElOaO) 
END 


C  95 
C  96 
C  97 

-  C--  9B 

C  99 
C  100 
C  101 
C  102 
C  103 

C  IDS 
C  106 
C  lOT 
C  lOi 
C  109 
_£_U0— 
C  111 
C  112 
C  113 
C  U4 
C  115 

C  111 

-  C  IIB 
C  119 
£  120 
C  121 

-£-122— 
£  123 

C  125 

C  127 
_X-_L2B_ 

-X— 130-  - 
£  13l« 
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OUTPUT 


SUBROUTINE  OUTPUT  J  I 

IHPLICIT  REALMS U-HtO-2) ^  -  _ _ _ _ 

COMMON  /OEN/  PEI » AMI «AME#DPDXf PREF <2) »PR (2) »P <2) »DEN» AMU?XU*XO*XP»  U  2 

2IMEN®IHEATfZ®TO#lNTe/V/U<EO0) 9F<i #200) #R <200 ) t RHO <200 ) fOM<200) *y<2  J  4 

30QJ/C/iC<2Qa)  ®AU<2O0)  tflU <200)  ^U <2MJ  9Ai24L2AlU-tB12t200)  »C(2»20aiyO-  U _ 5- 

4/YR<2aO) fUR<200) fRR<200) fHR<200) 9XM<200) ®PITOT<200) «TEMP<200)/E/OS  J  6 

5TARC30O)  iXRSOOO)  tRWRSOOO)  ®COSAi„<3O0)/l/N«NPUNP2»NPJjN£<UWPHiK£X  U  _  7 

§®K|NfKASE®KRAD/B/BETA®6AMA<2) tTAUI #TAUE» AUI <2) «AUE<2) •INDI <2)9lN0E  J  6 

JXEJ _ J _ 9— 


l/F/ITURBfIPRINTf TURBIN® TURSFSflDEUY 

IF  liNTOsNEa)  00  TO  i  1  _ _ _ _ _ _ J-  10 

00  14  Ia|tNP3 

.  YR  ill  ®00D#Q  _ _ _ _ _  _ 

ALaALPHA®iaO#0«0/3a41S92650*0  J  11 

16® It  CQMCI>-®l@l«NP3)  —  _ _ _ _ 

WRfTE<6®6)KRA0®I0IHEN®NEQ*KIX®KIN9lHEAT*ITURB®|0ELY®IPRINTfN 

XRSC300)8O®D^Q  _ 

CONTINUE  U  14 

If  CDFLQAT  < INT6®l)/0FU0AT  <IPRINT) ®NE®DFLOAT  t  C INTO-l )/lPRINT) JHETllRN _ 

DPDX@aDPDX^OAH/<@AM»l®D«0)  U  15 

IF  <KRAolE@rO)  @0  TO  2 

OS  sR  < I ) ®OSORT  <R  < 1 ) #®2«290^0®CSALFA®  <R<ll®Y iNP31»0.SQ«<l»CS 

IALFA®Y<NP3)^^2-PEI/<RH0<NP3)®U<NP3) ) ) ) 

.  If  <CSALFA®N£®0®OOtO)  OS  sOSTAR (IMTaX/CSALFA _ 

CONTINUE  J 

_nQ-3_X®UNRa _ U_ 


J  1? 
J  19 


21 

-21- 


TENP  <1 ) 1 1 f I ) -U  <  X ) #®2/2®O^0 
--  JUUI)  alKll  /OSQRT  t  tOAM»l  *0^0)  »T£MEtl4 1 

KS@aXM<l)®^2 

_ ixiici)  ax®i#am04  -Riioi4XiJi  i 

^COAM-lsO^O) ) 

_ l2--Xmtl4L®U^U4mA0|i  GO  TO  3 


J  24 
__U  25 
J  26 
._27 


P1T0T<1)*<  <OAM^I*O#0)#XiQ/2eD^O)«®<6AM/<SAM»l®D»O) ) 

PITOT  <  ll*PlTAT-m#44^AMa.®n^l/I2^a^0#mM«XSa»<IAM^UOm  liaU02 
^<®AM»l®0#0) ) 

_ COMTXMUE-1. _ _ _ _  - _ _ - _ _ _ 

S@2«OSQRT <2®D«0) 


PITOT < I) aPITOT <l) /PITOT <NP3) 

UR<I)aU<l)/U<NP3|..  - 

RR<I)«F<g®l) 

HR<I|sU<l)/SQ2 _ _ 

CGNTINUE 


J  29 
a  30 

U--11 
J  32 
-J — 33- 


U 

O- 


34 

-36- 


-J 

U 


3a 

39 


00  iS  Kt«3®N 

.  ._I^IAsTM£TA^.(T£MP.tMP3iy-T£MPiKU  tOR4KU-®  <  UO.#-0®UR-tfU.).X^ _ _ _ _ 

I  TeMPCNP3)/T£MP<KU^l)«UR<Kl.«U®<leO^O»UR<KL«n  )  )/2®0#0®<Y<Kl.*U® 

.  . —  -  -  ■  - -  - 


IS. 


THETAbTHETA^TEMP  <NP3) /temp  <3) ^UR  <3) ®  < I ^D^O^UR  <3) ) /2®O«0»V (3) ♦ 

a - ;£MiUi4gim£MRiNmAaR4iiie44^^  <  np  a  (np3  )  «y  < 

»NP1)) 

.  -..Cr*TAUlXWR4NPll-«^PHOiAIP114- _ 

aRITE  <6*8) 

UIUftaXAUIOfifiES^ _ 


J  40 
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OUTPUT 

1  USTAROOO)  «XL 

WRITE  (6«13)  AJI(n»OS  »P£1 «AME*0P0X6«DX»XRS (300 ) •USUP>PPO» 

1  CF. THETA 

WRITE  (6*9)  U  44 

DO  5  Isl*NP3*2  J  45 

WRITE  (6,U)  (Y(I)  .F(Itl)  .URd)  »HR(I)  .Hd)  •RHO(I)  tYRd)  •XMd)  tTEMP  J  46 

I  d)  .RRd)  »PITOTd)  )  J  47 

5  CONTINUE  J  48 

WRITE  (6*10)  (Y(NP3)  .Fd*NP3)  .UR(NP3)  *HR(NP3)  *R(NP3)  •RHO(NP3)  *YR(N  J  49 

1P3) *XM(NP3) .TEMP(NP3) •RR(NP3) .PIT0T(NP3) )  .  J  50 

WRITE  (6*9)  J  51 

RETURN  .J  52 

C  J  53 

6  FORMAT(////////»20X*»THE  INPUT  FLAGS  ARE» .///*25X# »KRAD  a  •*11*  J  54 

l//*25X* MDIMEN  «  • *  1 1  * // » 25X , ♦ NEQ  *  ♦ *  11 *//*25X» »KEX  »  t*!!*  J  55 

2//*25X»»KIN  s  NII.//.25X.  dHEAT  »  ♦  *  1 1  • // .  25X .  ♦  ITURB  »  »*Il* 

»//*25X*»I0ELY  «  »*I1* 

3//*25X**IPRINT  «  »*I3*//25X*»N  =  •*13) 

7  FORMAT (24H1THE  VALUES  OF  OMEGA  ARE/ (IPlOEl 1 *4) )  J  57 

8  F0RMATdMl*4X*»INTG**9X*'XU»*nX*»Z»*9X**RE0HS»*7X.«GAMMA»*aX*»PR»  J  58 

1 *9X*  »ZETA»  *ax*  »TURBI»*7X*  »TURBE»  *5X*»0ELTA»  »9X* 

1  ♦XL»/*6X*  »QW»*9X*»DS 

2TAR»*8X*»PEI»*9X*»AME»*8X*»DP/DX»*0X*»OX»*8X**  REX  ♦ * 7X ♦ » USUP » • 0X 
3*»P/P0»*9X*»CF»*8X*dHETA»/) 

9  FORMAT(6X**Y»*9X*»H/HO»*9X*»U/UE»*8X*»U/UM»*9X*«R**10X*»RHO»*10X*  J  62 

I  »YL  ♦  *9X*  »M*  *UX*  »T»  *10X*  *E/HO*  *7X*  ♦PITOT*/) 

10  F0RMATdPllE12,5*/)  J  64 

11  F0RMAT(1P11E12*5)  J  65 

12  FORMAT(I8*4X*IP10E12.4)  J  66 

13  FORMAT(  1P11E12.4*//) 

END  J  68» 
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SOURCE 


SUBROUTINE  SOURCE  (J»l9CS90S) 

IMPLICIT  REAL«8(A-H#0«Z) -  - - -  - -  - 

COMMON  /6EN/  RE  I « AMI « AM£ >OROX vPREE (2) »PR (2) •P <2) vOEN* AMU»XU«XD«XP »  A  3 

■~IXU7DX,  XSTEP  f  USALTATALPMAl  XR .  REORS  t  GAM  vZETAf  PPO  *TWTO“,  YST  ART  tUSUP » ID - A — W~ 

2IMEN? IHEAT*ZsTO»INTG/I/N»NPl »NPatNP3*NEQtNPH*KEX»KlN*KASE»KRAD/B/0  A  5 

3^Tff7GAMAT^)T^AUTf  rAUE.A'JI  (2)  » A'UEIZ) VINOITl)  UNDI(2) /V/U (2011)  .F(2.2 — S - S“ 

400) »R(200) 9RHO(200) •OM(200) »Y (200) /C/SC (200) fAU (200) 9BU(200) »CU(20  A  7 

50 )  V A  f 2 f  200 )  9 B  ( 29  200  ) » C  (21 200 ) /O/YRTIOir)  # URl^OO r»PR  f200yVHRT200T7XM  A - 0“ 

6(200) 9 PI  TOT (200) 9  TEMP ( 200 ) /E/DSTAR (300 ) 9XRS (300) vRWRS (300) tCOSAL (3  A  9 

TOO)  ■  ■  “  --  . -  . .  - 

IF  (J«6T#1)  GO  TO  1 

- C5*sc  ( IT » ( orr#  d  wr  i  ♦  i )  »u  ( i )  ( i ) )  /  (  om  ( i  ♦! )  -om  ( i ) ) - 

CS«CS»SC(I-l)»(U(I)»U(I)-U(l»l)«U(I»l) )/(OM(l)®OM(I-l) ) 
C$®(KD«0»UO^0/PR£F(J)r«CS/^(DMn>l)-OM(r-U)  - 

DSs0«O«0 

■  return 

1  CONTINUE 

- -  I  FTTaFr^ )  QQ'IXTT - 

REX»OSQRT(PPO«»( (l90«0-GAM)/6AM)-l,D*0)«PPO»»( (6AM^l,D^0)/(29D40<»6 
TAM)*  Z£TA«(GAM"T«0^y7GW)  ' 

REXsR£ORS»XO»REX 

XRS(300)fiREX  - 

CONs3a0O«0«U82O-09<»REX 

IF  (CON7tT72«TD«0)XON82'*TD«0 - 

3  CONTINUE 

■”  DELTAsOSTAROOO) 

ETAfiYd  )/DELTA 
IF{rTA*U£.0»2040) YLsYTTr 

IF (ETA9GT90,2D40.ANO«ETA«LE.O*4D*0) YLaY(I)/2,O*0*O,lO*0*DELTA 

rFn(FrA7iriTOT^lI^0  9ANO.ErA9LE.0«5D40)YL®0,3O40«yn)y(57l8D^0#DELTA - 

IF(ETA9GT,095D«0*ANOtETA*L£.0«60*0) YL«-0.1D*0«Y(I) ♦0.38D*0*DELTA 
IF(ETA«6T«0,6D*0»  AND«  Ef AVUE  ®  0  «  70*  0 ) YU=»0 • 20*  0»Y ( I )V0  .>40^0*OEUT A 
IF (ETA»6T,0»7D*09AND«ETA.LE.0»aO^0) YLa“0o4D«0*Y(I)»0.58D*0«OELTA 
IF  (£TA,6Ted,0D4.O)YLs  (-Y(T)*l9424D^O«D£LTA)  /2*40«0 
yR(I)aYL 

- IF  (l9tO,UOR,I.EQ.NP3)“CS*O.T5*0 - 

IF  (I.EQsKOR.UEQoNPS)  60  TO  2 
CAUL  VEFFdtl^iyEMU) 

CSs  (Ud^l)«-Ud-1)  )/(OMd^l)»OMd«l)  )»Rd)/P£I 
CS«CS«»0«320*0<»RHOd)»F(29l)  - 

CSbCS-EMU«F  (29l)/(RHOd)»U(I)  »YL»YL)»C0N 

2  CONTINUE 

EEEs®32D#0»P£I/(Ud)»Rd)»(  (Ud«l)®U(I-l)  )/(OM(Ul)-OMd-l)  )  ) ) 
‘EEE®EEE^CON^F { 2 9 1 ) / ( YL«YL»RHO d ) »U (I ) )  ■ 

DS!«CS/F(29l)-£EE 

RETURN 

END 
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READY 


SUBROUTINE  READY  E  I 

IMPLICIT  REAL»8(A®H#0®Z)  . .  . 

COMMON  /GEN/  PEI  *  AMI  * AME*DPDX»PREF (2) sPR (2) .P (2) #DEN» AMU«XU»XO*XPf  E  2 

lXL»DX»XST£P>CSALFAiALPHAtXR>REORStGAM»ZETA«PPOtTyLI{UlSIAfiTillSUP-aD _ A__4_ 

2IMEN9lHEAT»Z»TO*INT6/V/U(200) tF {2f 200) *R (200) fRHO(200)  ♦OM<200)  »Y (2  E  i 
300)/I/N9NPlfNP2fNP3fNEQfNPH9KEXfKIN9KASE#KRAO/B/ittAj6AMAmtTAUl4  5  . 

ATAUEtAUI (2) »AJE(2) tINOI (2) 9lN0E(2)  E  6 

CALL  DENSTY  _  E  I__ 

call  rad  (XU»R(U tCSALFA)  E  @ 

_C. _ Y  NEAR  THE  I  SOUMQARY _  _ _ _ E__^t_ 

GO  TO  (192*3) 9  KIN  E  10 

1  Y(2)®(KD«04.iETA)«OM(3)«»4.D«Q/t(3,D->0«RHQt2)tRHOl3)  )»(Ut2)^U(3)))  E  11 

GO  TO  4  E  12 

2  Y(2)®129D«0<»OM(3)/((3®D»0»RMO(2)«RHO(3))«(U(2)«U(3)«4,DtO»U(l)))  E  13 

60  TO  4  E  14 

J3 _ Jim  ®^5D^i)»0Mxaixtfii4aiii»u  id  j  _ e _ is_ 

4  Y(3)aY(2)«»25O«0«OM(3)®(l,D^0/(RHO(3)»U<3) )«29O«0/(RHO(3)»U(3)^RHO  E  1§ 

l(2)»U(2)))  E  17 

C  Y  »S  FOR  INTERMEDIATE  GRID  POINTS  E  IS 

DO  S  I®49NP1  E  19 

5  Y(I)aY(I“l)^.5D^O»(OM(I)®OM(I«l) )<»(l9D«0/(RHO(l)®U(I) )*l.O*0/(RHO(  E  20 

_ ll-D^^UXI®!) ) )  _ _  -  _ _  _ £  ZX 

C  Y  NEAR  THE  E  BOUNDARY  E  22 

Y  (NP2)  aY  (NPl  )♦  #250*0®  (OM (NP2)  -OM  (NPD  )  »  (l#D*0/  (RHO(NPD  ®U  (NPD  )  *2#  E  23 
lD*0/(RHO(NPl)®U(NPl)^RMO(NP2)®g(NP2) ) )  E  24 

GO  TO  (69790)9  KEX  E  25 

Y(NP3) aY(NP2) *(l,O*0*BETA)®(OM(NP2)»OM(NPl) )®4#D*0/( (RH0(NP1)*3#D*  E  26 

4(l®RH0iNP^jLX®lUtNP-L)*UXi4PaiD  -  - _ 

GO  TO  9 


6 

7 

0 

-9- 

C 


_ E-  27 


y  (NP3)!iY(NP2)*l2®0*0®  (0M(NP2)»0M(NPD)/(  IRHO(NPD+3#D*0®RHO(NP2) )® 

1 (U(NP2)*U(NPl)*4o0*0»U(NP3) ) ) 

GO  TO  9 

Y (NP3)«Y(NP2) ♦#5D*0® (OM(NP2)®OM(NPl) )/(RH0(NP3)«U(NP3) ) 

_l£_iCSALFA#Efl«iUD*G#nR^KRA£)#^£Q^O)  GO  tO_D _ _ 

XXX  IS  USED  TO  KILL  SHOT  IF  NECESSARY 

DO  10  Is2iNP3  _ 

XXX»R(l)»R(l)»2#O*0®V(I)®PEI»CSALFA 
IF  (XXX«LT#0#0D*0)  XOs2#O*0®XL 
IF  (XXX»LT#090D*0)  60  TO  14 

L»P£4 /  <R  ( I )  *OSQRT  (R  <  D  tft.(  I )  «#2»D*L0ftl-U-LtJ^EX®£SALFAXI — 


£ 

E 

E 

E 

E 

E_ 


26 

29 

30 

31 

32 
_33 


E  34 
E  35 
E  36 
E  37 
E  38 


c 

11 

12 

13 

CHANGED  SIGN  OF  2  IN  THE  DENOMINATOR 
-  GO  TO  13  . 

DO  12  |a29NP3 

Y(I)»P£I®Y(I)/R(D  . - . 

Y(2)®2#0*0®Y(2)«Y(3) 

Y  (NP?)  i«2*D4.0»¥  (NP2)  »¥  (NPD 

OF  ABOVE 

FOR  INTERNAL  FLOW  E 

f 

40 

41 

42 

43 

44 

45 

■  -  - 

E 

- . .  .  _ _ E 

E 

E 

C 

CALCULATION  OF  RADII 

E 

46 

_  . 

DO  14  I®29NP3  _  -  . 

_  _ 

-  .  - . . £ 

47 

IF  (KRAD#£Q#0)  R(I)®R(D 

E 

43 

_ _ 

IF  (KRAO#NE»0)  R  ( I )  aR « D  «Y  (D ®CSAtFA- 

,  ^  ... 

. . . . . E 

49 

C 

CHANGED  SIGN  IN  EXPRESSION  ABOVE  FOR 

internal 

FLOW  E 

50 

1 4 

CONTINUE 

E 

.M 

IF  (R(NP3) .LE#0#00*0)  XD»290*0»XL 

E 

52 

IF  (Y(NP3)^LT«0^00*0)  XD«2#D*0®XL 

_ _ _  _  E 

53 

RETURN 

£ 

54 

-END 

E 

55^ 
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VEFF 


SUBROUTINE  VEFF  (I9IRS9EMU) 

IMPLICIT  . 

COMMON  /6EN/  PEI  9  AMI  9 AHEtOPOXoPREF (t) tPRd) 9P(2) tOENt AMU9XU9XO9 XPi 
1  XL  9  OX  9  X  S  TEP  9  CSALF  A  9  ALPHA  9  XRj  REORS  t  SAN  @  lEI  Aj  PPOjJJiLias-TSIAaiAiiSUP^iJ- 
ilMENtIM£AT9l9T09lNT6/V/U(aOO) tF (ataOO) 9R(200) 9RHO(200) 9OM(a00) 9Y(a 
300)/I/N9NPl9NPa9NP39NEQ9NPH9KEX9KIN9KAS£9KRAO 
T®Fa9i>-u(X)»®a/a90«o 
TT^FClf  IPI)»U(IPn»«^a/2tO«0 
T®(T<^TT)/a9D«0 

EMU^T^^ZETA/  (REORS/DSQRT  la-aD^Qli _ 

IF  <I.EO®UORtKEQ®NP3)  RETURN 
EEsO®3D^O®RHO(l)^F(afl)^PEI 

EE^EE/ ( RMO ( I ) «U < I) ^R  U ) ^ (U a « 1) “U < I “ 1 n / (ON ( S ^ i ) -ON ( I- U ) ) 

EMU®EMU«E£  _ _ _ _ _ 

RETURN 

END 


P  I 

p  a 

A-  -  4.  ... 

P  4 

P  S 

P  4 

P  7 

P  8 


P  10 
M — U^-- 


VISCO 


FUNCTION  VISCO  U) 

_  IMPLICIT  REALMS  (A-M 9  0-^ .  _ 

COMMON  /SEN/  PEI tAHI 9AME9DP0X9PREF (2) 9PR<2) 9P(a>  9OEN9 AMU9XU9XO9XP9 

XALAOXjA&I£RA£SALFA9ALPHAtXRj.R£ORSi6AN92£TA9i>PQtTWTQ9YSTART9USUPf  ID 

2IMEN9lHEATfl9TO9lNT6/V/U<a00) 9F  Utaoo ) 9R <200 ) tRHO (200 ) 9OM (200) 9 Y (2 
300) /I/N I NPl.i.Ng2AtjP 3tNEQtNPHtKEXfKIN 6 KA&EaKRAD  . 

TaFdf  I)-U(I>®#2/290«0 

IF  (l»EQ®l®0Rtl#E@9NP3)  RETURN 

_£i:«(UaD.*-j)»RHQ  (1)»F  LIaJUl®PU . . . . 

EEaEE/(RHO(l)®U(l)#R(I)i^(U(l«U-U(I-U)/(OM(i4l)-ON(I-l))  ) 
W!SCQ»V!SCQ»EE _ 

RETURN 

_EMO_ 


Q  I 

Q  2 
A  .4 
Q  4 
0  5 
Q  6 
Q  7 


FBC 


SUBROUTINE  FiC  (X9U9 INO9 AUFS) 

IMPLICIT  REALm(A-HjD«>U  _ : _ 

COMMON  /SEN/  PEI i AHI 9 AHE9DPDX9PREF (2) 9PR(i) 9p(2) tDENt AMU9XU9XDtXP« 

_ 1AL» QX.XSTEP.CSALFAtALPHA.XR.REQRStQAMtZETAtPPOtTMTQ.YSTART.USUPtID 


IND«2 

.,AafSa,0.aM.t.fi _ 

CONTINUE 

IN0«1 


H 


aiMENf IHEATfZtTOfINTS/V/UdOO) fF(2f200)#R(a00) fRHO(20O) 9OM(20O) 9Y(2 

300)  .  .... _ 

TW  IS  PRESCRI8EQ  IF  IHEAT  a  1  —  QO0T  IS  PREiCRISEO  IF  NOT  1  H 

-  iNO«i  .  _ _  ..: _ _ ^ _ fei 

AJFSsTWTO  h 

ClHEAT^EQa)  SO  TQ  I _ ^ _ tt 


M 

Jl_ 


2 

_4_ 

4 


S 


9 

-UL 


H  n 


RETURN 

END 


H  12 

-M-_X3j 
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COEFF 


SUBROUTINE  COEFF 

-  IMPLICIT-  BEAL _  _ _ _ _ 

COMMON  /OEN/  PEI t AMI » AMEt DPOXtPREF <2) tPR<2) fP(2) sOEN# AMUiXUfXOsXP* 
lXL«QXfXSTEP»  CS  ALF  A  »  ALPHA.  XR  «  REQRS  a  QAM  .  ZET  Aa^PPQ  (JMm.L¥SIARTjdiSUP-j^ 
ilHENf IMEATfZfTOf INT@/I/NfNP|fNP2fNP3fNEQ»NPMiKEXfKIN»KASEfKRA0/B/i 
aEIAaOAMAiZJ  ^TAliIj  tAUEAAUimAAJ£{2)  9INOI 12U INOE  (1X2^^ZU120Q1  jF.(2a.2 
400) 9RC200) fRHO(aOO) fOMtZOO) 9Yli00)/C/§C(200) »AU<200) fBU<200) 9CU(20 
§0)  @4i2t2O0)  iB(2tl00)aC(2r2©fll _ _ _ _ _ _ 


DIMENSION  Q1 (2QQ) t  62(200) t  63(2C 

13C200) 

CALCULATIDR  OF  _SMAU^jC_1S _ 

DO  I  I«2fNPl 

mm  9  5Dt  ft®LlRU  #  U»RUU _ 

RH®@iO^O®tRHOU*l)^RHO<l)  ) 

-UMa^CliLfttlUCItU  tUClU _ 

CALL  VEFF  Utl^ltENU) 

SC ( I ) bRA^RA»RH»UM#EMM/ (PElOPEl ) 
THE  CONVECTION  TERM 

SA®Bill^-AMl/F£I- _ _ 

SBs(R|NP3)^AME«R(U«AMI)/PEI 


00  4  I®3fNPl 

_  _ 

Pi®925D^0/OX 

F3®P2/0MD .  _  _ 

Pla«OM<l#I)®OHa)  )®P3 

_pj8  tmm®0Rii®iiJLtP3 _ 

P2i»390«0®P2 

QbSA/OMO _ _ 

R2»®Si«9iiO^0 

_ : _ _ _ 

RI®-(OM(l#l)«3»0«O«OMa))®R3 


B2JMMPZ±m _ 

63a)»P3«@«R3 

ggplf  ULLf  UjsPi 
THE  DIFFUSION  TERM 


__D _ 

D  6 
J 

D  S 
_D _ i-_ 

0  10 

0  12 

D  14 

_D _ IS_ 

D  13 

-_D _ ll_ 

D  IS 

0  20  " 
1 


0  22 
__D _ ia-_ 

D  24 

_Q _ ES_ 

0  23 

_Q _ 21— 

•0  2i 

0  30 

0  '32 


0  34 

0  33 

_D _ 31— 

0  38 


IF  (NEOsEOtl)  @0  TO  3 

HO  2-J»UMPM _ : _ 

C  ( Jf  1 )  U9 1 ♦  U  ®P2®F  Uf  I )  ®P3»F  C  Jf  I-»U 

f'Ai  i  emiors’  r .  i  .r'e .n  1  t .  f  t  1 


C  Uf  I )  •»€  Uf  I )  #CS®F  (Uf  I )  ®0 I ) 

-MJ-OJjmjJJjmEE-LU _ 

BUtD^euiD/PREFCU) 

--C0HIXNUE _ _ 

SOURCE  TERM  FOR  VELOCITY  EQUATION 

Cl  I  f  I  Rni»n»#n 


S2  a ) 8P2®il a ) / <  RMO  < I ) ®U  C I ) ) 
SUI)®Pl®SUI)/IRH0a^l)»U(l^l) ) 


D  43 

_ D _ 41_ 

0  4a 
_ D _ 49^- 

0  50 

h  sii 


0  St 

P  54 
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COEFF 


Si  (I)^Sl(l)/U(I«U  0  S6 

S2(l)pS2(I)/Ua)  D _ 5Z 

S3{nBS3(I)/U(I-l)  0  SS 

CQNXIMJi^ _ _ D _ Sa. 

COEFFICIENTS  IN  THE  FINAL  FORM  D  60 

DO  5  liiJsNPi  D _ tX- 

RL®l«0^0/(G2(I) ♦AU<I) «BU(I)®S2(n )  D  62 

AU(1)«{AU(XJ  ♦SKD^iai  0 _ 6a 

BU(n»(BU(I)^S3<I)-e3(I)  )»RU  0  64 

-3^ _ jCLLLii  0 _ 6a 

IF  (NEQtEQeU  GO  TO  7  0  66 

DO  6  D  67 

DO  6  Ia3fNPi  0  68 

RL^l9Q»0/(6am>AXJtII»AiU4X)«D{4aXJ _ D _ 6a 

A(J9l)8(A(JfI)®61(l))«RL  0  70 

_ B  {  y  G3JLLU-aRL _ Q 7X 

6  C{Jf  l)^CU9l)®RL  D  72 

7  CALL  SLIP  D _ Za 

RETURN  0  74 

END  _D _ 7a 


DENSTV 


SUBROUTINE  OENSTY  F  | 
IMPLICIT  REAL«a(A-H#0»2)  _  _  _  _ 

COMMON  /GEN/  PEI » AMI  tAMEtOPDXf  PREF  (2)  »PR  (2)  »P  (2)  tDEW,  AMU»XU*XO»'xp*  F  2 
IXLtDXfXSTEPtCSALFAiALPHAiXRtREQBStGAMtZETAtPPOtTWTQ.YSTART.USUPtlD  A  4 
2IMENf IHEATf Zf TOf INTG/V/U(200) * F ( 2® 200 ) * R (200 ) » RHO { 200) »OM{200) *7(2  F  4 


300)/l/NfNP|  fNP2»NP3fNEQ?NPHfKEXfKIN»KAS£iKSAD. _ _ _ F  ^ 

TNP3sF(lfNP3)-U(NP3)»«2/29D«0  F  6 

RHONP3«TNP3««^(l90*0/(GAM™l90*0)  )  F  7 

DO  1  I«l9NP3  F  i 

T^Fdf  I)-U(I)»<»2/2®D^0  F  q 

1  RMO(I)»RHONP3^TNP3/T  F  10 

RETURN  _ F  11  . 

END  F  |2-» 
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ENTRN 


SUBROUTINE  ENTRN  0  % 

_ IMPLICIT  REAL»B(A-HfQ-2)  _ 

COMMON  /OEN/  PEIf AMIeAMEfOPOXfPREFCa) iPR(a) jP<a> sDENiAMUfXUtXD^XP*  6  2 

_ LXL « DXAj^I£Pjd:LSAJLi^^  xa.9  REORS  9  QAM « ZETA  «.P.P.QjaJimUJtST.JLRT.j  USUP « in _ i 4^ 

ilMEN9lHEAT9l9T09lNT@/V/U<aOO) 9 F < 3# 300 ) iR (300 ) t RHO (200 ) »OM (200 ) e Y (2  @  4 

_ 3QQ)/l/N9NPl9NP2iNP39NEQ9NPH9KEX«KIN8KASEtKRAD _ (2 _ S _ 

1  60  TO  (2f396) 9  KEX  66 

_ ^2 _ RETURN- _ a ,  7— 

3  CONTINUE  6  8 

_ C _ TM£L_FQLLQMIN6  AME  IS  FOR  LAMINAR  FLOtrf _ 6 _ 9_ 

IF  (INTOsNEti)  60  TO  5  -  6  10 

_ 00  4l_  I*itNP3 _ 6 _ U _ 

IF  (OM(l) 96T«0»90«0)  N9®I  6  11 

_ IF  (QM(I)  ,6T«0.9Qt0)  ImNRJ _ 6 _ L3 _ 

4  CONTINUE  6  14 

_ OOLMaOM(N9)»OM.iN3»U _ 6  15 _ 

DOM9*QN(N9)®09904-0  6  16 

U9®U(N9)-D0M9/00M»(U(N9)-U(N9®U )  6  IS 

_ TERMBu  (U  CN9)  "U  (N9-| )  )  /DOM _ 6 _ _ 

R9»R(N9)«D0M9/00M»(R(N9)®R(N9®1) )  6  10 

_ RH09«RH0  CN9>!!U20M9/DOM»  (RHO  (N9L«RM0  tNOalU _ 6 _ 2X— 

VIS9^VISCO(N9)«OOM9/OOM»(VISCO(N9)«-VISC0(N9-1) >  6  21 

_ OUUPgR  tN9l  »RiN9)  #RHO  lN91«UJlN9jJ»VXSO(2  tN9) _ 6_  ^ _ 

CUUP«CUUP^R9^R9®RH09«»U9»VIS9  6  14 

^ _ CUUEaCUUPi-t  2^*0  »P£  1 1 _ _ 6 _ 15 _ 

CguM^R<N9®l)#R(N9®i)«RH0(N9-«l)®U(N9®U«»VISC0(N9-l)  6  16 

_ £yUMgCUUMiR9»iR9tRH09t.U94  V I S9 _ 6  21 _ 

CUUMwCUUM/(2®0«O^PEI )  6  28 

_ _ - _ Q5al^JljU}»CUUPy  (DOM4DOM9) _ _ 6 _ 1® _ 

66®2tD^O®CUUM/(OOM«(099D'8>0«»OM(N9'»U  ) )  6  3© 

_ _ T£aMAaflS*-CUm9)HU9X-66»(U9"U(N9-U  ) _ 6 _ 31 _ 

TERMA®TERMA/TERHB  6  31 

_ U®MlJ3gUauMJQSaflaLia^a4  0 1  ^OSQRT  (  I  ^0»©LttPPQ»t  |  <Q4MyUD40I/6AMU _ _ 6 _ 33 _ 

C  USUP  IS  REAOoIN  IN  SEOIN  IT  SUPPRESSES  THE  i«L9  6  34 

— _ _ TERMCw  (U9WIG«U9I/DX*DgOX/  {RH09»U9) _ Q _ IS _ 

T£RMC®TERNC^PEI/TERNi  6  36 

_ AMEMlERMAmXERMC _ _ _ 6 _ 37 _ 

AHE®AME-®09lD^O4R(l)4AMI  6  38 

- AM£«AME/(0^9D*0*R(NP3) ) - S _ 39 _ 

RETURN  6  40 

— 6 - - 0 — 44 

RETURN  6  42 

- EN0-  -fl _ 63a_ 
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MASS 


SUBROUTINE  MASS  (XUfXDeAH)  I  i 

IMPLICIT  REALMS (A-H«0-Z)  _ _ 

c  applicable  to  an  impermeable-wall  situation  I  a 

.  AM^O.D^O  !  1 

RETURN  I  4 

END  ...  _ I 


PRE 


4 

5 

6 


SUBROUTINE  PRE  (X®OPOMK)  K  1 

IMPLICIT  REALMS CA-HfiO-I)  _ _ _ _ _ 

COMMON  /GEN/  PEI t AMI ©AMEePPOXfPREF (g) iPR <a) f P (a j t OEN» AMU? XU»XD»XPf  K  2 

_  ilLt  DX.t^IEF  @  C  SALE  A®  ALPHA  j  XB  aa£QaSteAM®ZETAaPPQtIMTQ®YSTAaijdLiS4i£iLlD _ k A- 

21MEN»IHEAT®Z@TO@INTG/V/U(2O0) ®F (a®aO0) f R (100) 9RHO(aOO) »OM(aO0) ®Y(i  K  4 

30 Q )  /1/Nf  NP 1  ®  NP2  f  NP3  ®  NEQ  9 NPH 9 KEX  t  K IhU KASE t KRAP _ IL.. 

DIMENSION  XX(300)®  POP(30O)  K  A 

IF  C1NT6®NE,0)  go  to  1  .  _  _ _ K_..  1... 

read  (593)  LMAX  K  8 

_  REAP— ISaAI-  IXX  CL)  *PQPiL)  aL®!  tUHAXi _ X _ S_ 

WRITE  (Gf5)  K  ■  10 

WRITE  (6@6|  (XX  (LI  ®PDP  (t)  eLaL^lfLHAX)  _ -K.__ai- 

CONTINUE  K  la 

Lai  _ _ K.  -03-. 

CONTINUE  ■  K  14 

_ _ K — IS- 

IF  (XX(L)tLT.X)  60  TO  2  K  1§ 

DPDXa(POP(L)-POP(L'^l))/(XX(L)®XX(L®l))  _  _ K_a7_ 

PP0»POP(L-I)^OPOX»(X-XX(L«»1))  K'  18 

DPDXX«OPOX  _  _ _  _  K.  19 

OPDXXaOPOX«(6AM»KO^O)/GAN  K  20 

-  -  AEJHRN. _ ^ _ _ _ _ _ K _ 


F0RMAT(13) 

FORMAT(2E12@0) 

FORMATClHl@/saX@  »X®  al2X-®  «P/P0»9/J 
F0RMATClP2E|4t5il6) 


K 

K 

K 

__X- 

K 

_X_ 


22 

.£3 

24 

-25- 

28 
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RAO 

SUBROUTINE  RAD  (X*R1 *CALPHA)  U  1 

implicit  REAL»0(A-H*O”Z) 

COMMON  /GEN/  PEI .AMI *AME»DPDX.PREF (2) .PR (2) .P (2) •DEN« AMU*XU*XD*XP»  L  2 

IXL  .DX.XSTEP.CSALFAf  ALPHA  ♦XR»REORS»eAMiZETA«PPQiTWTQxlSIMJ^USUP»iD 
2IMEN.IHEAT9Z*TO.INTG/V/U(200) .F(2.200) *R(aOO) ♦RHO(200) »OM(200) .Y(2  L  A 

300)/I/N*NPI»NP2»NP3»NEO»NPH.KEX9KIN»KASE.KRAO  L  5 

applicable  to  nozzles  with  CONSTANT  LONGITUDINAL  RADIUS  OF  L  6 

CURVATURE  OF  THE  CONVERGING  SECTION-  CONSTANT  WALL  HALF  ANGLE  OF  L  ? 

DIVERGING  SECTION-ANO  WALL  SLOPES  MATCHED  DOWNSTREAM  OF  THE  THROAT  L  8 

IF  (INTG.NE.O)  GO  TO  1  _ _ L _ 9.^ 

PI2s3«14l59265D*0/2.D«0  L  10 

ALPHAaALPHA«PI2/90.D*0  L  11 

COSALFaDCOS(ALPHA)  L  12 

SINALF®DSIN(ALPHA)  L  13 

ZWiGsXR* ( I .D^O^SINALF)  L  U 

XWIGkXR«(PI2^ALPHA)  _ 1 _ 15 

RWlGal.D«O^XR»(UO^O-COSALF)  L  16 

1  CONTINUE  L  17 

IF  (X,GE«XWIG)  GO  TO  2  L  18 

R(1)b1.D*0^XR»(1.O^0-OSIN(X/XR) )  L  19 

CSALFAsOSIN(X/XR)  L  20 

CALPHA^CSALFA  _  _ _  _ ^_L _ 2i_ 

R1=R(1)  L  22 

ZZsXR»(l,D*0-OCOS(X/XR) )  L  23 

GO  TO  3  L  24 

2  CONTINUE  L  25 

R(1)®(X-XWIG)«SINALF*RWI6  L  26 

CSALFA®C0SALF  _ _  _  _ _L_  27 

CALPHAsCSALFA  L  28 

RlsR(l)  L  29 

ZZ«<X-XWIG)»CSALFA»ZWIG  L  30 

3  .CONTINUE  L  31 

!F(!D|MEN.£Q«0)RlaUO«0 

ZZsZZ-XR  L  32 

ZsZZ  L  33 

RETURN  L  34 

END  L  35- 
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SLIP 


1 

a 


SUBROUTINE  SLIP  M  I 

IMPLICIT  REAL«»8(A-H#0»Z) 

COMMON  /6EN/  PEI « AMI t AMEsDPDX iPREF ( 2) sPR (i) »P ( 2) *DEN* AMU ^XUsXD.XPf  M  2 
iXLiJ)XaXST£P»CSALFA#ALPHA#XR#R£ORS»SAM«ZETAtPPOiTWI(UYSTART»USUPtID  A  4 
2IMENf IHEATf Z#T0@INTS/I/N»NPUNP2»NP3*NEQ*NPH*KEX»KIN9KASE»KRAD/V/U  M  4 
31200)  tF(Ef200)  fRiaOO)  8RHO(200)  iOM(20Q)  »¥  Uftai^B/BOAtOAMAlE)  iTAUla.  M^  5 
4TAUE9AUI (2) 9AJE(2) f INDI (2) iINOEia)  M  6 

COMMON  fU  AK»ALM6/C/SC(200) «AU(200) tBUiaOQ) »CUi2Q01fA12aa00) iB{2a  M  7 
1200) «C(2f200)  MS 

_ _SL1P  COtf  F 1 C I E N T S  N£ Afi  Ji4£  l.BOUNQARY  FOR  -VELOCITY  EQUATION _ At. .  .3 . 

CU(2)^09D*0  M  10 

CUlNPD^OeD^O  _ _ _ _ _  M  n. 

GO  TO  119283) 9  KIN  M  12 

BU(a)aO®0«0  _ _ _  _ _ _  M  13 

AUia)sUO«O/lKO^O«2,O^O»0ETA)  M  14 

.  60  TO  A _  .  .  ..  _ _ _ M..  .15 

M  16 


SQs84,D«0®U(l)«Ull)®l2.D^0«Ull)«^U13)«9,D«0»U(3)*U13) 

iU12)®88D«0®12,D«0®Ull)^U13))/t2*  D  ♦  0  ®1U1J  ♦X.DtMO  U1  .♦OSORT  iSfl)  I  . . 

AU(2) sl®O«0«§U12) 

60  TO  5  _  -  -  _ 

BU12) s0#0«0 

-  CALL  VEFF  l2*a^£MU)  _ _ 

AK 1  sKD^O/OX-OPOX/  (RMO  1 1 )  ®U  ( I )  ®U  11 )  ) 

AK2«®UU)®AKl#DP0X/lRH01i)®U(l) )  _ _ _ _ _  _ 

AU«RHOU)®U11)®®25040®1Y12)*Y(3)  )®®2/EMU 

IF  1KRAD,£Q®0)  60  TO  4  _ _ _ _  .  - 

AU12)®2«0*0/(2®0^0«AU®AKl) 

CU12)  5®s§0*0®AJ®AKli»AUlEJ  _  . _ _ _ 

GO  TO  5 

CU12)  sUD«0/12®D<90#3®0^0®AJ®AK1)  _ _ _ 

AU<2)®CU{2) ®12tD^0-AJ®AKl) 

CU12)a«»CU(2)®4®O»0®AJ®AK2  . . 

SLIP  COEFFICIENTS  NEAR  THE  E  BOUNDARY  FOR  VELOCITY  EQUATION 

60  TO  16f798)  9  KEA  _ _ _ _ _ 

AU1NP2)®0@0*0 

BUlNP2)«i*O#O/U9D»O^2®D»Q®0ETA)  _ _ _ _ _ — _ _ _ 

60  TO  9 

SOa84®O^0®UlNP3)®U(NP3)-12«D*0«UlNP3)  »U  tXlPU^9L,a*A4UXMPlJ-«UtNPll 
AUlNP2)a8«0«0® 12®0«0®U1NP3) ^U INPl ) ) / 12»0®0®U (NP3) ♦T.D^O^U INPl ) ♦DSQ 
_*AT<SQ)  ) 


M 

M 

M 


17 

18 
If 


M  20 
M  22 


M 

M 

M 

M 


23 

24 

25 

26 


M-— 27-- 
M  28 
M  29 
M  30 
.  M.  31 
M  32 
— M _ 33_ 


M 

_M_ 


34 
3a 
M  36 
M  J37 
M  38 


0U<NP2)®UO«O-AU1NP2) 

60  TO  9  _  _ _ _ 

AU1NP2)®0®0'90 

CALL  VEFF  INPl  ®NP2*£MU4 _ _ _ _ _ 

BKl»l®O«0/OX®OPOX/<RHO(NP3)®U(NP3)®UlNP3) ) 
_Bii2M^  1MP3  )®BKl40PDX/(RHQ<NP31»iJlNP3)  > _ 


M  39 


M  41 
M . 43- 


M  43 
_M _ 4A_ 


9 

_C- 


0U®RHO  1NP3) ®U 1NP3) »*25D®0® 12®0«0®Y 1NP3) -Y INPl ) ®Y 1NP2) ) »«2/EMU  M  45 

CUlNP2)®l«0«0/(2*D#0^i«i)AA«aO®aiaJ _  Jl_  Jk6_ 

BUlNP2)®CUlNP2)®12®0^0-iU®BKl)  M  47 

CUlNP.21®®CUlNP2)®4«0.-9j)tflJ®BK2  _ _ _ _  ..  48 

IF  1NEQ®£Q®1)  RETURN  M  49 

SLIP  COEFFICIENTS  NEAR  THE  I  BOUNDARY  FOR  OTHER  EQUATIONS _ M _ 51L 


00  20  U®l9NPM 
ClJiL2)®0.9DtQ 


M  51 

M_.  52._ 


C1J9NP2)«00O^O 
_  60  TO  CIO i  1254.3)  I  KIN 


M  53 
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SLIP 

10  CALL  FBC  (XD»  J9 INOK  J)  »Qn  M  55 

IF  (INOl  (J)  .EQ.D  60  TO  U  H  56 

AvJI(J)3QI  M  57 

AUt2)®l.D«0  _  _ _ _ I4_  Sfi 

B(Jf2)s0«D4>0  M  59 

C{j9  2)s8,D^0»(UD*0«2*D*0OBETA)*PREF(J)  »AJI  ( J) /  UK«AK»fl£Tl«  ( M  60 

1BETA)®(1,D«O^BETA)4»{3.D^O<»RHO(2)*RHO(3))«U{3))  M  61 

60  TO  15  M  62 

U  F(J*l)sQI  M  63 

A(Jt2)»{l.D»0»B£TA«6AMA(Jl)/(UQtO»B£TA^GAMACJJJ _ — _ 

B<j92)i«l9D^0»A(J»2)  M  65 

60  TO  15  M  66 

12  A(J?2)®(U(2) ♦U(3)®8tD«0»U(l) )/(5,D^0O(U(2)*U(3) )^8,D^0«U(1) )  M  67 

GFs(KD«0“PREF(J)  )/(l.D^OtPREF(J)  )  M  68 

A(J92)«(A(J92) ♦6F)/(l,D^0*A{J.2)«6F)  M  69 

B(vJ92)*8a#DtO"A(J*21  ^  _  - _ _ M  74 

GO  TO  15  M  71 

13  B(J»2)»0*D40  M  72 

CALL  SOURCE  (J»l*CS»0S)  M  73 

AK1®I.D40/DX«DS  M  74 

AK2s»AKl«F(J»l)«CS  M  75 

_AUF«AJ»PR£F{J)  M  76 

IF  (KRAD,EQ,0)  GO  TO  14  M  77 

A(j92)s2.D«0/(29D*0^AJF«»AKl)  M  78 

C(J*2)a'»®5D*0«AJF»AK2»A(J«2)  M  79 

60  TO  15  M  80 

14  C(J#2)®1«D«0/(2,0«0*3,D^O»AUF»AK1)  M  81 

A(U»24aC(j921«(2*0«O^AJE»AKl) _ _  _ -M-_a2- 

C(Jf2)®“C(J»2)«49O*0«AUF«AK2  M  83 

C  SLIP  COEFFICIENTS  NEAR  THE  E  BOUNDARY  FOR  4TMEH  £aUATIONS  M  84 

15  GO  TO  <16» 18»19) »  KEX  M  85 

16  CALL  FBC  ( XO  t  J*  INOE  ( J)  9QE)  .  _ _  N  86 

IF  (INDE(J)  ,EQa)  60  TO  17  M  87 

AUEtJlsQE  N 

8(j9NP2)«l9D«0  M  89 

A(j9NP2)sO«DoO  M  90 

C(J9NP2)«»8,D«0<»(1.D*0*2»0*0<»BETA)  »PREF(J)»AUE(U)/(AK^AK«0ETA*(UD  M  91 

1^04BETA)«(l.D4Q*BETA)»(RHQ(NPl)*3,DtO»RMQ{NP211*U(NPlU _  M  92 

GO  TO  20  M  93 

17-—  X(J»NP31sQ£ _ _ _ _ J1__J94- 

B(U»NP2)a(l*D«0^BETA-GAMA(J) )/(1,0+0«8ETA«GAMA(J) )  M  95 

A(U9NP2)3UD^O’-0(J*NP2.)  N  96 

60  TO  20  M  97 

18  S(j9NP2)8(U<NP2)^U<NPn»8,D*0<»U(NP3)l/(5,D'>0»(-UlNP2)-*yLNP4rn*8*D*a  M  98 

««^U(NP3)  ) 

- OEia  1.^0^  a®RR£E_UU^  iy4^^0^IUP-R££-U44_ _ M-  99 — 

B(U»NP2)s(8(J»NP2) ♦GF)/(1.D*0+B(J9NP2)#GF)  M  100 

A(J9NP2)8UD*0«»B(U»NP2)  -  -  -  MlOl 

GO  TO  20  M  102 

19  A(UfNP2)s0.O#0  -  - . M  103 

CALL  SOURCE  ( J»NP3»CS*DS)  M  104 

-  .  .  8Kl8l,D*0/0X-DS  _ _ _ _ _ -H-445- 

BK2®«0K1»F(J»NP3) »CS  M  106 

BUFsBU<^PREF(J)  M  107 

C(JfNP2)®UO^O/(2.D^O^3.D«O»BUF»0Kl)  M  108 

B(JfNP2)«C(U»NP2)«(2-D*0«ajF»BKll  _ _ _ _ -  .  M  109 
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SLIP 


C(J9NPE> a-C<JfNP2)«4.D^O®BJF»0Ka  M  HO 

aO  CONTINUE  _ M_1J4 _ 

RETURN  N  112 

END  -H  Xia^ 


SOLVE 


SUBROUTINE  solve  ( At B»C»F *NP3)  N  1 

IMPLICIT  REAL«8(A<“H«0»Z) 

C  THIS  SOLVES  EQUATIONS  OF  THE  FORM  N  2 

C  Ftl)  s  A(I)«F(l«i)  ♦  B(I)«F(1-1)  ♦  C(I)  _ _ _ ..„X 

C  FOR  l®2»NP2  N  4 

DIMENSION  A(NP3),  B(NP3)»  C(NP3)»  F(NP3)  N  5 

NP2aNP3-l  N  6 

B(2)sB(2)»F(l)^C(2)  N  7 

DO  1  I33«NP2  N  0 

T®1«0^0/(1.D*0”B(I)«A(I»1) )  N  .  .9 

A(I)«A(I)®T  N  10 

1  B(I)®{B{I)<^B(I»1>  ♦C(I)  )«T  N  11 

DO  2  I329NP2  N  12 

JsNP2-I^2  N  13 

2  F (J)®A(J)^F(J«1)*S(J)  H  U 

RETURN  N-  15 

END  N  16» 


WALL 


SUBROUTINE  WALL  R  I 

IMPLICIT  REAL«>0(A«H»O-Z)  _  _ 

COMMON  /6EN/  PEI # AMI « AMEfOPDXt PREF (2) t PR (2) ♦P (2) »DEN» AMU*XU9XD*XP»  R  2 

2IMEN#IH£AT*29TO#INTG/V/U(200) »F(2’l200) 9R(200) »RHO (200) *OM (200) *7 (2  R  4 
300)/I/N#NPl»NP2fNP3.NEQ»NPH»KEXfKIN#KASE»KRAD/fl^SE^tAr6AMI121^IAUlt  _8  5 

4TAUEfAUl (2) f AUE(2) •INDI (2) »INOE(2)  R  6 

C  CALCULATION  OF  BETA  FOR  THE  I  BOUNDARY  ...  _  . _ _  R  7 

1  YI®«5D^0«»(Y(2)*Y(3)  )  R  0 

UI®»5D«0<»(U(2)  ♦UO)  J  S _ 9 _ 

RHiS925D^0»{3»D»0®RHO(2)  oRHOO)  )  R  10 

RE«RH»UI*YI/V1SC0(1)  _  _ _ _ _ _  R  11 

FPsOPDX»Yl/(RH<»UI®UI)  R  12 

AM®AMI/(RH®UI)  _  .  _  R  13 

C  FOR  LAMINAR  FLOW  AND  AM=0  (NEED  DIFFERENT  EXPRESSION  IF  F»0)  R  14 

SpI«D«Q/R£-FP/2^D^0  --  - _ _  _ .8 _ IS _ 

B£TAsRE«(S*FP«AM)  R  16 

TAUI«S«»RM®UI«UI  R  17  . 

IF  (NEQtEQel)  RETURN  R  10 

C  CALCULATION  OF  GAMA  »S  FOR  THE  I  BOUNDARY  ..  .  .  R  19 

DO  2  U=1»NPH  R  20 

C  _ FOR  LAMINAR.  FU3W  .AhUl  LIMIT INa.  ZEflIL.AM _ B _ 21 _ 

SF®1«D^0/(PR(J)»RE)  R  22 

6AMA(U)»RE»PR(J)<»(SF^AM)  .  ..  R  23 

IF  (INOI (U) *EQ*l)  AUI (J)sSF»RH«»UI«(2#0*0»F(J»l)-F(J»2)-F(U»3) )»»50  R  24 
»40  _  -  _  _  _ 

2  CONTINUE  R  25 

_ RLTURM  _ B _ 26 _ 

END  R  27® 
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The  following  is  an  example  input  to  the  boundary-layer  computer  code. 


See  subroutines  BEGIN  and  PRE  for  input  formats. 

00321004150067 

2,880E 

05  ,667  .85  1,4  ,0  ,0 

642. 

,9999 

.1 

1.0  .1  .85  .01  .000001 

022 

0,0 

.5283 

32,1 

.4390 

64.2 

.3609 

96,3 

.3012 

128,4 

.2533 

160,5 

,2151 

192.6 

.1850 

224.7 

,1612 

256,8 

.1470 

266,9 

,1370 

321. 

.1278 

353.1 

,1220 

385.2 

.1128 

417.3 

.1077 

449,4 

.1027 

481,5 

•  1003 

513,6 

.0965 

545,7 

,0950 

577.8 

,0943 

609,9  .0935 

642, _ ,0935 

9000.  .0935 
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MAIN 


C  THIS  IS  A  DATA  REDUCTION  PROGRAM 

C  INPUT  IS  PTl?  TW/TO»  XMEf  GAHMAf  REYNOLDS  NUMBER  PER  FOOTt  AND 
C  PT2  DISTRIBUTION 

C  UNITS  OF  PTl  AND  PT2  MUST  BE  CONSISTENT 
C  The  INPUT  XME  IS  FOR  A  FIRST  GUESS  ONLY 

C  IF  THE  INPUT  TW/TO  IS  ZERO  THEN  THE  WALL  IS  TAKEN  AS  ADI0ATIC 
C  OUTPUT  IS  VELOCITY?  TEMPERATURE?  BOUNOARY^L A YER  THICKNESS?  BOUNDARY- 
C  LAYER  DISPLACEMENT  THICKNESS?  BOUNDARY-LAYER  NOMEMTUM  THICKNESS? 

C  shape  FACTOR?  MACH  NUMBER?  TBAR?  AND  SKIN  FRICTION 

C  HEAT  flux  can  ALSO  BE  DETERMINED  (SEE  THE  CORRESPONDING  LAB  REPORT) 
DIMENSION  Y(IOO) 9PT2(100) ?PT2P1 ( 1 0 0 ) ? XM ( 1 0 0 ) 9TOTOE(100) 

DIMENSION  U{100)?T(100) ? TBAR (100) ?RHO(100) ?RHOU(100) 

DIMENSION  CF  (100)  tCCCdOO)  ?CC(100) 

_ Ct-  (1)  a»001  _ _ _ _ _ 

RtA0(5?2)NUM 

_ DO  10  LL«1?NUM _ 

READ(5?l)  PTl ?TWT0?XME?6?R 
REA0(5?2)NE 

READ (5? 3) (Y{N) ?PT2(N) ?N=1?NE) 

C  RtAD(5?l)  PTl ?TWT0?XME?G?R 

C  READ(5?2)NE 

C  REA0(5?3) (Y(N) ?PT2(N) ?Nsl ?NE) 

1  format (8E10,0) 

2  FURMAT{I3) 

3  FORMAT (2£l0e0) 

_ GPlsG^U _ 

GM 1 =G- 1 s 

Ds  (GPl/2,  )  <»<»  (6/GMl ) 

IF (XME,LE«1«) XMEsSQRT {2./GMl«( (PT1/PT2(1) ) »» (GMl/G) « ) ) 

_ XE  =  XM£^XME _ 

IF  (XMEeGTd.)  CALL  XEDGE  ( PT2  ( NE  )  ?  PT 1  ?  G  ?  XME  ?  XE ) 

PI=PT1/(1,*GM1»XE/2,)»»(G/GM1) _ 

PT2(l)sPl 

_ X^XE _ _ 

XM(NE)»SQRT(XE) 

DO  4  N=1?NE 

4  PT2P1 (N) =PT2 (N) /PI 
J  =  NE-1 

DO  5  N^1?J 

_ IsNE-N _ 

IF  (PT2P1  (I)  .LE«D)  X  =  2./GM1#  (PT2P1  ( I )  <>«^  (GMl/6) -1  ®  ) 

_ IF  (PT2P1 (I) .GT.D)  CALL  XDIST (PT2P1 (I) ?G?X) _ 

XM(I)sSQRT(X)  . 

5  _ CONTINUE _ _ 

XME=XM(NE) 

A=GMl«XMt<*XME 
CALL  PR(XME?G,A?P) 

_ HW  =  TWTO»  (  K»A/2«  ) _ 

Cisl* ♦A/2«-HW 

_ IF  (TWTO»EQ«0. )  C1  =  0, _ _ 

D=SQRT  (2»»A»HW>?C1*»C1) 

_  IF (TWTO.£0>0»)DsSQRT(A»2>»(l.*A/2»)  ) _ 

U(1 ) =0, 

_ DO  7  K=1 ?NE _ 

C=(XM(K) /XM(NE)  )»<»2 

CALL  VT(A?HW?C1?D?P?C?T(K) ?U(K) ) 


91 


AEDC-TR-76-62 


main 

c"i^L  TEMP  ( A  9  H  W  f  cVf’Q7(^U"('K )  « DHD  V  t  T  ( K )  ) 

IF(K.EQ.l)T(K)sH»y _ 

1bAR(K)=(T(K) ♦A«U(K)»U(K)/2.-HW)/(l,«A/a9»HW) 

_ KhO(K)=l,/T (K) _ 

RnOU (K) sRHO (K) (K) 

_ iXC  ( K)  =1  .-RHOU(K) _ _ _ 

CC  (K)  sRHOU  (K)  <  1  .«U  (K)  ) 

_ _ lU T OE  ( K1 = J T  ( K)  ♦A<*U(K)»U(K)/2.)/<l,»A/2.) _ 

IF  (K,NE,1)  CF(K)=CF (K-1) 

_ IF (K.NE.l)CALL  SKiNCRtCP  (K)  t6f XMEf Y (K) tU(K)  ) _ 

T  (K^n  =T  (K) 

7_  continue _ _ 

PT?Ptl»pf2(NE)/PTl 

CALL  SUM(YfCCC»NE»VINTl) _ 

CALL  SUM(Y»CC*NE,VINT2) 

_ CALL  CFI (XME»R»V1NT2»SQCF) _ 

CALL  PLUS  (SOCF,  Y»RHO*R»NE •CCC*CC* VINTl * VINT2* A*HWf Cl » 0*P) 

_CFF=SQCf- «SQCP_ _ 

SMAPE=VlNTl/VINf2 

15  M=l9NE _ _ 

IF (U(M) «GT. 0*99)60  TO  16 

15  CONTINUE _ 

16  DtLTAsY(M»Ut(0.99-U(M«l)  )»(Y<M)-Y(M»1)  )/(U(M)®U(M-l)  ) 

_ JF  J  TWIO  «  EG  ._0  .J  W«  lTE(6tll)  P»6»PT2PTI  tRtOELTAt  VINTl  tVINT2tSNAPEtCFF 

11  F0RMAT<1H1  *39Xf  »"thIS  DATA  REDUCTION  SHOT  IS  FOR  AN  ADIABATIC  WALL* 

_ 1  »//  »40Xt  <RECQ VERY  FACTOR  s  0.8800  *  t//t40Xt  »PRANDTL  NUMBER  g  U 

2  F6.4»//»40X» *6AMMA  *  » »F6.4,//»40X» »PT2E/PT1  *  • » 1PE10.4*//* 

_ 3  40Xf  »FR£E-STREAM  UNIT  REYNOLDS  NO.  PER  FOOT  «  «  » lPE10.4»//t _ 

4  40Xf ’DELTA (IN)  =•*  IPEl 1 .4f//*40Xf ’DELTA  STAR(IN)  «•»  lPE11.4»//» 

■  5  40X» ’THETA(IN)  s’t  IPEl 1 .4t//>40Xt »SHAPE  FACTOR  «»#  lPE11.4»//» 

6  40X9 ’CF  (BASED  ON  HE-THETA  EQ  OF  R.  &  M«  REP.  NO.  3712)  *»* 

_ 7  1PE11.49///) _ _ 

IF (TWTO.NE.O.) WRITE(6»12)P9G9PT2PTl9R»OELTAfVINTl»VINT29SHAPE9CFF 

12  FORMAT(lHl939X9 ’THIS  DATA  REDUCTION  SHOT  IS  FOR  A  NON^ADI ABAT IC  WA 
ILL’ 9//940Xf ’RECOVERY  FACTOR  =  0,0800 ’ 9//94OX9 ’PRANDTU  NUMBER  *  ’9 

.2  F6.49//94OX9 ’6AHMA  s  »  9F6.49//»40X9 ’PT2E/PT1  ■  ’ 9 lPE10.4t//f 

3  40Xf ’FREE-STREAM  UNIT  REYNOLDS  NO.  PER  FOOT  «  » # IPEIO ®4f //* 

4  40X9  ’DELTAdN)  IPEl  1 .49//94OX9  ’DELTA  STAR(IN)  «’  9  1PE11.4»//. 

5  40X9 ’THETAdN)  =’9  IPE 1 1 . 4  9 // 9  40X  9  » SHAPE  FACTOR  ®’9  lPEll.4»//9 
_ 6  4QX>»CF  (BASED  QN_RE-TH£TA  EQ  QF  R.  I,  M.  REP.  NO.  3712)  »»9 

7  iPEll.49///) 

_ )tllIXE16  t0J _ 

8  FORMAT(  /97X9’Y’9llX#’PT2’98X9’TO/TOE’9lOX9«M’9lOX9’U/UE’9  9X» 

_ 1  ’T/ TE<  9  9X9 ’TBAR’ 97X9 ’RHO/RHOE’ 92X9 ’RH0»U/ (RH0E»UE) ’94X9’CF’9/) 

WRITE (696) (Y(N) 9PT2(N) 9T0T0E(N) t XM (N) t U (N) 9T (N) 9TBAR (N) 9RHO (N) 9 

_ IRHOU(N) 9CF(N) .Nel.NE) _ 

6  FORMAT  dP10El3,5) 

WRITE(69l8) _ 

10  FORMAT (///) 

10  CONTINUE _ 

STOP 
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Xe.DG£ 

SUBR"OU TWE~XFDGfc  (WtZ» GTyT^T) 

_ XaY^Y _ ^ _ 

6P1=G4.1,  ' 

Gm1=G-1« 

Csw/Z 

1  F«(GP1<»X/(GM1#X«2^)  )«^»(G/GM1)»(GPI/(29^G«X-GM1)  )<»<»(  K/GMl) -C 

DFDX"  (F4-C)  *6/GM1»  (  2«/ (X^  (GMHX-9-2#  )  ) -2®/ (2«»G»X«GMn  ) 

XNEW=X“F/DFDX 

IF  (ABSnx-XN£W)/XNEW)  sLTsOsOOODGO  TO  2 

XsXNEW 

GO  TO  1 

2  X=XNEw 
HtTURN 
END 


xorsT 

SUBROUTINE  XUIST (B*G»X) 

GP13G4-1  9 
GMlsG-K 

1  Fs(6Pl^X/2@)^«{G/GMl)<^(GPl/(2»»G^X-GMl)  )»^(1«/GM1)-B 
DFDX»(F4B)*G/XM2««X'“17)/{2®<^6«^X-GMI1 
XNEWcX-F/DFOX 

IF  <ABS(  (X-XN£W)/XNEW)  «LT®0«000UGO  TO  2 

XsXNEW 

GO  TO  1 

2  X=XNEW 

return 

END 


PR 


SUBROUTINE  PR(XM,GtA?P) 

C  CALCULATES  PRANDTL  NUMBER  FOR  A  RECOVERY  FACTOR  OF  0,08 
Rs,88 

AL=5»»7®/2« 

ALlsAL^U 

AL2sAL^29 

BEs4»^AL 

As  (6-1  •  )  «XM«^XM 

DsSQRT  (2»<»A<»  (U^A/2*  )  ) 

CALL  TINYF ( A«HW»0» tDf Oe 9FO) 

CALL  TINYF (AfHWf 0# fDf 1 e 9FI ) 
hONE=A«BE/ ( AL1»AL2) ❖A/Zs^Fl-FO 
PsK»a^{R-U)/(2»«H0NE) 

RETURN 

END 
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skIn 

BR 0“UT r NE " SKT N  HI? TO  7^1^ 9 XM »Y*U) 

Cl=b,7b 

CH=S«10  .  ■ 

0MEGA=b768 

REsR<>Y/12« 

ASSORT {CF/2« ) 

Z=T#-^«88»  (  (GAM«l  «)  /2«^ 

_ ^^SQRT(Z^  _ 

S1GMA=  {  ■(  (GAM»r,  )72« )  7  (1  (  (  (‘GAM»1  «  )  /2.  )  ^XM>«2  )  ) 

SQSIG=SQHT (SIGMA) 

1  GAsCl<»AlOG10(HE»X<»Z<*<»("*“2e^(®5^0MEGA)  )  )  *C2 
G  =  GA-  ( 1  «  /  (SQSIG^X^Z)  )  AHS I N  ( SOS  I  Gi»U ) 

DgT)X=  (C1/(2;3^*X)  )7  ( 1  a/  (SQS1G<»  (X*»2)  ^T)  )  <*ARS  IN  (  SQS  IGT^U  ) 
XNEW  =  X”(5/0GL)X 

IF  ( A B S  {  (  X N E W « X  ) 7 X Nl; w  )  « LT.  0  «  0 yo“iT  G 0  T 0^  2 

X=XNEw 
GU  TO  1 
^7»XNEW_ 

CF=27^X^“'»2 
heturn 
t  NO 


SUM 

S  U  B  R  0  U  T  i  N  E  S  U  M  (  xTY  f  N  E  9  \flW)~ 

C  X  IS  THE  INDEPENDENT  VARIABLE 

C  Y  IS“ THE' DEPENDENT  VARIABLE 

dimension  X  ( 100) f Y  ( 100) 

V1NT=0* 

NN  =  NE«r 

'  DU  1  Isi ,NN 

1 _ ViNT  =  VlNT^(X(I^l)»X{I)  )^{Y(I4-1)«Y(I)  )/29 

RETURN 

END 
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PIUS 

SUBROUtlNt  PLUS(SQCF*  Y » RHO ♦  H * NE ♦  CCC ♦  CC .WNT l*VlNT2*A»HWtCl»U»P)  ^ 
DIMENSION  YdOO)  ♦RHO(IOO)  ♦CCCdOO)  *  CC  ( 1 00 )  •  YP  ( 1  00  )  » UP  ( 1 00 ) 
DIMENSION  CCCP( 100) »CCP( 100) *  TP n 00) 

C  THIS  SUBROUTINE  USES  THE  WOLL  MODEL  TO  DETERMINE  Y  AND  U  NEAR  THE 

C  V»ALL  TO  IMPROVE  THE  CALCULATIONS  OF  THETA  AND  DELTA  STAR 

Yt  =  138,/RH0{1)«SQRT(2*)/SQCF<»12./R 
DO  1  I=2»NE 
1F(Y(I) .GT,YE)GO  TO  2 

1  CONTINUE 

2  N=1 

DY=Y (N) -YE 

CALL  SUM (Y9CCC*N,VINT1W) 

CALL  SUM  (Y,CC9N,VINT2W) 

VINT1=VINT1-VINT1W 

VINT2=VINT2-VINT2W 

YPLUS=0. 

YP(1)=0. 

UP(1)=0. 

CCP(D=0. 

cccp(n=i. 

DO  3  J=2*70 
YPLUS=YPLUS+2. 

UPLUS  =  ATAN(8.<»TAN(0.09#YPLUS/8.)  )/0.09 

YP  ( J)  =  YPLUS/RHO  ( 1 )  <»**1  •268<>SQRT  (2.  )  /SQCF#12,/R 

UP  (J)=UPLUS<»SQCF/SQRT(2.<»RHO(l)  ) 

call  TEMP(A,HW9C1»0*P*UP(J) 9DHDV*TP<J) ) 

CCCP( J)=1.-UP( J)/TP{J) 

3  CCP(J)=UP<J)/TP(J)»(U-UP(J)  ) _ 

UvJ®  7  0 

_ CALL  SUM(YP9CCCP» JU»VINT1P)  . _ _ 

CALL  SUM(YP»CCP* JJ»VINT2P) 

V  INTI =V INT !♦ VINT  1 P+ (CCCP(JJ) ♦CCC(N) )/2.*DY 
ViNT2®VINT2>VINT2P+  (CCP(JJ)  *CC{N) 
return 
End 
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CFI -  - 

SUHROuriNt  Ch I <XM£ VlNT2f SOCF ) 

This  SUriWOUTlNt  solves  f- OH  Cf-  USIinG_  tOUAT  ION  (22)  OF  THE  telINTER 
AI\D  GAUOtr  HfcPOHT  M«  NO,  3712  OECtMHER  1970  ' 

_  1^0=1 

t-C  =  SQHT  (!“•♦,  2^XMt<»XME)  - - 

C  =  f'0<»R/12,oVlNT2 
0=0,3894 
,  n=4.632 
XNFw  =  0,002<^*0,b 
bb  X  =  X!\fcw 

F  =  1,0-"H<*X  ~  '■ 

TtXP=txP (0.537/X) 

XNE  w  =  X~  (  (C-i*F*^TtXP“C)  /  (G<>F  '^TEXP<*  ( ”0  •  537/X»«2 )  “H<>TEXP)  ) 
b  =  ABS  (  (XnFw-X)  /X!mEW) 

IF (S,LT, 0,0001 ) 00  ro  100 
bO  TO  bb 

loo  ■  CONTINUE  ' 

5oCF=XNEw 

hltukn 

End 


The  following  is  an  example  input  to  the  data  reduction  computer  code. 
See  MAIN  for  input  formats. 


001 

17,68 

0,  ,8  1.4  3819000 

021 

,0 

13,53 

,00bb 

14,71 

,034 

14,94 

,07  7 

15,18 

,11 

15.28 

,177 

15.75 

.206 

15.9 

,2492 

16. 

,2778 

16,11 

,3206 _ _ 

16.16 

,349b 

16,37 

.3924 

16.41 

,421 

16,54 

•  4641 

16,62  . 

,507 

16, 6J 

•^357 

J6,69  ^  .  . _ 

,5786 

16.77 

jl6216 

16.92 

,65 

17,07 

.6932 

17.16 

1,4 

17,60 
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A 

a 

C 

c  ^ 

Cf 


Ch 

Cp 

Cy 

E 

e 

e 

Fc 

Fs 

f 

H 

H 


H 

h 


h' 


h 

ho 

hi 


NOMENCLATURE 

(7  -  1)M^2 
Constant  in  Eq.  (6) 

Defined  by  Eq.  (13) 

Constant  in  Eq.  (48),  equal  to  4 

Local  skin-friction  coefficient,  liXjlpJJj') 

Stanton  number,  -  qw /p«,Uoo(haw  -  h^) 

Specific  heat  at  constant  pressure 
Specific  heat  at  constant  volume 
e/U,-2 

Turbulent  kinetic  energy  per  unit  mass,  (<u2>  +  <v^>  +  <w2>)/2 
e/Ho.. 

(1  +  0.2  M„2)1/2 

1  +  0.056  M_2 

Function  defined  by  Eqs.  (60)  and  (67) 

Specific  total  enthalpy,  also  used  for  shape  factor,  6*/0 
H/h^ 

Specific  enthalpy 

Fluctuating  local  static  enthalpy 

h/h^ 

Zeroth-order  h 
First-order  h 
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k  Molecular  thermal  conductivity 

kt  Eddy  thermal  conductivity 

L  Dissipation  length,  L  =  L/s 

£  Prandtl  mixing  length 

Integral  scale  length  of  turbulence 
M  Mach  number 

m  "u  =  77!/^,  used  herein  as  m  =  7 

m  Mass  flux 

Pr  Prandtl  number,  /x  Cp/k 

Pr^  Mixed  Prandtl  number,  Cp  (/x  +  fXt)/(k  +  kt) 

Prt  Turbulent  Prandtl  number,  Cp  jXt/kt 

p  Mean  pressure 

p  (t  - 

Po  Impact  pressure 

q  Heat-transfer  rate 

Red,c  Reynolds  number  based  on  pipe  diameter  and  velocity  at  pipe  center 

Reynolds  number  based  on  pipe  diameter  and  mean  velocity 

Reg  Reynolds  number  based  on  boundary-layer  thickness  and  free-stream 

velocity 

Reo,s  Reynolds  number,  Po  s/po 

Re_^^x  Reynolds  number,  U^x/p^ 

r  Recovery  factor,  r  =  (Taw  -  TJ/(To,^  -  TJ 

f  Defined  by  Eq.  (5) 
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S 

T 

T 

t 

U 

U 

Ur 

u 

u 

u+ 

-<uv> 

V 

V 

V 

w 

X 

X 

y 

y 

y^ 

z 

a 

P 


Reynolds-analog  factor,  s  =  Cf/2ch  (also  used  for  characteristic  length) 
Mean  static  temperature 

(To  -  Tw)/(To,^  -  T^) 

Time 

Mean  velocity  in  x-direction 

U/(Ho,Ji/2 

Friction  velocity  (r^/pw)^^^ 

Fluctuating  velocity  in  x-direction 

u/u„ 

U/(Tw/Pw)1'2 

Reynolds  stress,  a  bar  over  this  quantity  indicates  normalization  with  respect 
to  Ho 

u  ^oo 

Mean  velocity  in  y-direction 
V/(Ho,Ji/2 

Fluctuating  velocity  in  y-direction 
Fluctuating  velocity  in  z-direction 
Coordinate  along  body  surface 
x/s 

Coordinate  normal  to  body  surface 
y/s 

PwUtV/Mw 

Coordinate  completing  right-hand  system  with  x  and  y 

A  constant  equal  to  0.18  in  Section  3.0  and  equal  to  (5/2)m  in  Section 
4.0 

c(5/2)m 
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7  Cp/cv 

A  Defined  by  Eq.  (62)  for  constant  temperature  wall  and  by  Eq.  (69)  for 

an  adiabatic  wall 

8  Boundary-layer  thickness  defined  as  the  value  of  y  where  U/U^  =  0.99 

5*  Boundary-layer  displacement  thickness 

e  A  constant  equal  to  zero  or  1/8  in.  Section  3.0  and  equal  to  (1  -  Prm) 

in  Section  4.0 

^  Defined  by  Eq.  (61)  for  constant  temperature  wall  and  by  Eq.  (68)  for 

an  adiabatic  wall;  also  exponent  of  T  in  power  law  relation  for  viscosity 

V  y/s 

6  Local  wall  angle 

6  Boundary-layer  momentum  thickness 

fjt  Molecular  viscosity 

M  m/Mo 

jUt  Eddy  viscosity 

p  Density 

P  P/Po 

Ue  Schmidt  number  for  turbulent  kinetic  energy 

r  Total  shear  stress 

Tt  Turbulent  shear  stress 

0  Stream  function 

CO  Transformed  stream  function 

V  nip 
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SUBSCRIPTS 


aw 

E 

I 

o 

w 


Adiabatic  wall 

Outer  edge  of  boundary  layer 
Inner  edge  of  boundary  layer 
Local  total  conditions,  except  ho 
Wall  conditions 
Free-stream  conditions 


SPECIAL  NOTATION  AND  CLARIFICATION 

<  >  Indicates  time  average,  e.g.  <uv>  = 

2t 

The  variables  denoted  in  the  nomenclature 
understood  to  be  time  averaged. 


(uv)  dt' 

-t 

as  mean  quantities  are 
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